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TWO -PHASE  PRESSURE  DROP 
A  LITERATURE  SURVEY  AND  CORRELATION  ANALYSIS 

Raymond  D.    Caughron 

September  1,    1966 

INTRODUCTION 

Any  combination  -  solid-liquid,    solid-gas,    or  liquid-gas  -  in  motion  con- 
stitutes two-phase  flow.     Recently,    considerable  effort  has  been  expended  on 
liquid -gas  flow,    since  this  type  of  mixture  is  vital  to  the  analysis  and  design  of 
refrigeration  cycles,   power  plants,    compressors,    condensers,    evaporators, 
boilers,    and  more  recently,    to  nuclear  reactor  research  and  space  rocket  technol- 
ogy.    In  these  investigations,    two-phase  pressure  drop  calculations  are  significantly 
important.     Unfortunately,    these  calculations  can  be  quite  difficult  and  complex  for 
the  following  reasons:    (1)  external  and  internal  bounding  surfaces  are  continuously 
interacting;  (2)  forces  and  thermal  interactions  at  the  interfaces  affect  changes 
in  the  fields  of  flow  velocities,    pressures,   temperatures,    and  thermal  and  diffusion 
fluxes;  and  (3)  twice  as  many  flow  and  property  variables  in  comparison  with 
single -phase  flow  theoretically  complicates  even  the  simplest  physical  flow  models. 

Because  simple  approximations  have  recently  proven  inadequate  for  most 
practical  applications,    emphasis  has  been  placed  on  methods  utilizing  physical 
models  correlated  with  experimental  data.     These  models  attempt  to  describe  the 
phenomena  without  becoming  so  mathematically  complex  as  to    render  the  method 
of  computation  used  impracticable.     Unfortunately,    the  methods  presently  avail- 
able still  are  too  complicated  and  inaccurate  for  consistent  application  in  standard 
engineering  equipment  design. 


This  report  summarizes  the  presenl  methods  of  two-phase  pressure 
drop  prediction,    gives  enough  essential  background  material  to  provide  a 
useful  physical  picture,  and,  hopefully,  provides  enough  details  to  stimulate 
additional  study  and  interest  over  the  field.     Additional  references  and 
excellent  reviews  have  been  published  bylsbin  et  al.  (28);  Jens  and  Leppert 
(30);  Tong  (48);  Bennett  (9);  Scott  (4:-!);  Anderson  and  Russell  (3);  and 
Gresham  et  al.    (20). 


LITERATURE  SURVEY 

A  true  scientific  basis  for  analyzing  experimental  data  on  single -phase 
fluid  flow  did  not  exist  until  1883,   when  Osborne  Reynolds  developed  from 

dimensional  analysis  the  equations  for  dynamic  similarity     and  dimensional 

2 

homogeneity.       Earlier  experimenters,   beginning  with  Couplet  in  1732, 

attempted  to  develop  empirical  equations  to  fit  their  experimental  data.   These 
equations  were,    however,  limited  in  that  they  applied  only  for  particular  pipe 
sizes  within  given  velocity  ranges  and  only  for  specific  fluids.     Following 
Reynolds,   in  1914,   Stanton  and  Pannell  (47)  investigated  single-phase  fluid 
flow  through  smooth  pipes,    applying  a  dimensional  analysis  to  the  fluid  friction 
problems  appearing  in  engineering  practice.     Their  experiments  on  compress- 
ible and  incompressible  flow  confirmed  the  feasibility  of  using  a  dimensional 
approach. 

In  1934,    serious  interest  was  taken  in  applying  the  theory  of  similarity 
to  the  motion  of  liquid -gas  mixtures.     The  following  excerpt  from  a  paper  by 
Schmidt  (42)  published  in  that  year  defined  the  situation  at  that  time: 
"The  motion  of  liquids  by  means  of  rising  steam  or  gas 
bubbles,    such  as  takes  place  in  steam-boiler  tubes  and  evaporator 
apparatus,   is  so  complicated  that  its  analytical  calculation  is 
impossible.   Evan  an  attempt  to  establish  the  necessary  differential 
equations  meets  with  insuperable  difficulties.     On  the  other  hand, 
resort  to  the  point  of  view  of  similarity  offers  the  possibility  of 
reducing  the  variables  to  a  certain  minimum  number  of  dimensionless 
magnitudes.     The  solution  of  the  problem  is  then  found  in  the  form 


Similarity  analysis  involves  dimensionless  numbers  with  selected  independent 


variables  used  to  describe  a  system  i.  e.,    Reynolds,  Weber,  and  Euler  numbers. 

2 

Dimensional  homogeneity  —  an  equation  is  dimensional 

the  dimensions  of  each  term  of  the  equation  are  the  same. 


4 
of  a  relation  between  these  magnitudes,    which  usually  have  to  be 
determined  experimentally." 

By  1944,    it  was  generally  realized  that  because  of  the  large  number  of 
variables  involved,    a  single  universal  equation  would  be  inadequate  to  correlate 
pressure  losses  over  the  entire  range  of  two-component  flow.     Martinelli 
et  al.    (35)  at  the  University  of  California  were  the  first  to  investigate  the  prob- 
lem,   using  Schmidt's  suggested  similarity  and  experimental  correlation  analysis. 
Their  studies  involved  isothermal  pressure  drops  for  two-phase  two-component 
flow,    without  mass  exchange  and  in  a  horizontal  pipe.     Eliminating  the  use  of 
flow  patterns,    they  obtained  a  fairly  satisfactory  method  of  solution  by  relating 
the  actual  pressure  drop  (APTp)  to  its  single-phase  counterpart  (AP   ,  AP  ) 
flowing  through  the  same  pipe.     These  studies  were  based  on  turbulent  flow  ex- 
perimental data  only.     In  1946,    Martinelli  et  al.    (37)  extended  their  studies  to 

the  viscous  region.     Later  Lockhart  and  Martinelli  (33)  further  extended  the 

3 
method  to  include  all  flow  regimes.       Various  investigators  have  since  attempted 

to  apply  the  equations  developed  by  Martinelli  and  his  co-workers  to  systems 

4 
with  and  without  mass  transfer  and  for  known  flow  patterns     in  vertical  and 

horizontal  pipes.     In  1948,    Martinelli  and  Nelson  (36)  adapted  the  two-component 

correlation  of  Lockhart  and  Martinelli  to  single -component  systems  by  a  simple 

correction  to  account  for  changes  in  the  axial  component  of  momentum.     The 

two-component  systems  of  Lockhart  and  Martinelli  did  not  involve  appreciable 

mass  transfer  between  the  phases,    whereas  the  single -component  systems 

of  Martinelli  and  Nelson  did.       Martinelli  and  Nelson  originally  did  not  possess 

sufficient   experimental   data   to   verify  their   approach.         Martinelli' s 


3 
Hereafter  referred  to  as  the  Lockhart -Martinelli  correlation.     See  page  43 

for  this  method  and  example  problem. 
4 
Methods  of  flow  classification  are  discussed  beginning  on  page  8. 


studies  provided  the  first  acceptable  methods  to  predict  two-phase,   two- 
component  or  one-component  flow  with  or  without  mass  exchange. 

In  1951,   Abou-Sabe   (1)   at  the  University  of  California   suggested  that 
the  Martinelli  correlation   (33,  36)  should  be  modified  to  compensate  for  the  effects 
of  flow  pattern  and  roughness  existing  at  the  gas-liquid  interface.     Baker  (7) 
later  developed  equations  to  calculate  the  pressure  losses  of  gas-oil 
mixtures  flowing  in  pipelines.     His  equations  essentially  modified  the  Lockhart- 
Martinelli  correlations  for  annular,    slug,    plug,    and  bubble  flow.     Prior  to  this, 
in  1946,   Armand(4),  and  Armand  and  Tretchev  (5)  in  Russia,  had  investigated  the 
flow  of  air-water  mixtures  in  pipes.     They  developed  three  equations  for 
correlating  pressure  losses  for  annular  mist,   annular  flow,    and  stratified 
bubble  flow.     Later,    studies  by  Gazley  and  Bergelin  (17),   Jenkins  (29),    and 
Holden  (23)  at  the  University  of  Delaware  showed  that  the  Lockhart-Martinelli 
correlation  was  not  applicable  to  stratified  flow. 

The  best  discussion  of  the  Homogeneous  method  in  which  a  two-phase 

friction  factor  (analogous    to  the  Fanning    single- phase  friction  factor)  is  used 

5 
has  been  presented  in  a  paper  by  Owens  (40). 

Harvey  and  Faust  (22)    later  combined  the  Homogeneous  and  Martinelli 
approaches  by  selecting  assumptions  from  both.     Rogers  (41)  then  utilized 
this  mixed-model  method  to  predict  the  behavior  of  two-phase  hydrogen  flows 
in  well  insulated  transfer  systems.     A  large  number  of  computations  were 
required  to  obtain  numerical    results  from  the  equations  of  Harvey  and  Faust; 
Rogers  used  a  high-speed  computer. 

Dukler  et  al.    (14)  from  his  "Comparison  of  Existing  Correlations" 


5 
See  page  39  for  the  Homogeneous  method. 

c 

Mixed  models  are  discussed  on  page   53. 


6 

7 
states  that  over  25  correlations  now  exist.     Dukler's  own  correlation     involves 

a  two-phase  Reynolds  number,    the  Koo  equation  (31),    as  well  as  other 
significant  terms  to  compensate  for  elevation,    including  an  iterative  procedure 
to  evaluate  the  acceleration  correction  factor.     This  method  predicts  the  two- 
phase  pressure  drop  independent  of  flow  pattern.     Additional  methods  by  Dukler 
utilizing  similarity  analysis  involved  the  following  hypothesis  (1)   assuming 
homogeneous  flow,    applicable  to  the  mist  or  spray-flow  region,    and  (2) 
calculating  a  two-phase  function  factor  and  Reynold's  number  in  terms  of  liquid 
properties.     These  criteria  are  applicable  especially  to  flow  patterns  with  high 
liquid-to-gas  ratios  such  as  plug  and  bubble  flow. 

o 

Chenoweth  and  Martin  (11)     have  introduced  an  empirical  correlation  for 

horizontal  flow  with  gas  and  liquid  superficial  Reynolds  numbers  greater  than 

2000.     This  correlation  used  larger  pipe  diameters  and  greater  pressures  than 

did  the  correlation  of  Lockhart  and  Martinelli  (33).     It  empirically  correlated 

the  ratio  (APTD/AP  ,)  of  the  two-phase  pressure  drop  and  pressure  drop  of  the 
IF         £* 

mixture  if  only  the  liquid  were  present,   to  the  liquid -volume  fraction  of  the 

feed  (R).     A  parameter  that  is  a  function  of  single -phase  friction  factors  obtained 

from  a  Moody  chart  appears  in  the  correlation. 

9 
Hughmark  (2  5)     empirically  correlated  the  ratio  of  the  two-phase  friction 

factor  (f^T-./f    .)  as  a  function  of  the  liquid- volume  fraction  in  the  feed  stream, 
iP    4* 

using  a  two-phase    Reynolds  number  as  a  parameter.         The  single -phase  friction 
factor  is  then  obtained  from  the  usual  single -phase  correlations,    using  average 
properties  for  the  gas -liquid  mixture  in  the  pipe.     Hughmark1  s  correlation 


7 
Dukler's  method  and  example  problem  may  be  found  on  page  33. 

o 

See  page  30  for  Chenoweth  and  Martin  correlation. 
9 
See  page  31  Hughmark's  correlation. 

ReTP  -  R<M<  +  R  u  (Eq.  25). 

srs.        g^g 


can  be  used  either  for  horizontal  or  vertical  flow.       Levy's  method 
proposed  in  1963  (32)  was  based  on  a  mixing-length  model.     This  model  treated 
the  two-phase  system  as  a  continuous  medium,   and  applied  to  it  methods  and 
assumptions  widely  accepted  in  single-phase  turbulent  flow.     In  1962,    Govier 
and  Omer  (18)  summarized  the  state  of  knowledge  concerning  the  calculation 
of  pressure  drop  and  other  related  quantities  in  two-phase  flow: 

"The  principal  flow  patterns  are  understood  in  a  qualitative  way  and  the 
effect  on  flow  pattern  of  the  major  variables,    mainly  the  mass  velocities  of 
the  phases,   is  recognized.     In  nearly  all  cases,   the  real  significance  of  the 
viscosities  of  the  phases,   the  possible  separate  roles  of  the  densities  of  the 
phases,    and  the  influence  of  the  diameter  of  the  pipe  is  not  known. 

"Extensive  experimental  data  confirm  the  fact  that  the  pressure  drop 
for  a  two-phase  system  is  influenced  by  the  flow  pattern,    and  indicate    the 
need  either  for  separate  pressure  drop  correlations  for  each  flow  pattern, 
or  for  the  incorporation  into  a  master  pressure  drop  correlation  of  these 
variables  which  in  fact  define  the  flow  pattern. 

"A  large  number  of  correlations  are  available  for  the  prediction  of  two- 
phase  pressure  drop  in  terms  of  the  major  variables  affecting  it.     Each  of 
these  has  its  own  shortcomings,   being  restricted  to  certain  flow  patterns, 
to  gases  of  specified  densities,   to  certain  diameters  of  pipe,   or  the  like. 
There  is  as  yet  no  generally  satisfactory  universal  correlation. 

Excellent  charts,   tables,    reviews,    and  comparisons  have  been  published 
in  the  literature,    namely  by  Bennett  (9),    Tong  (48),    Scott  (43),    Isbin   et  al. 
(27),    and  Gresham  et  al.    (20).     Experimental  data  (containing  over  9000 
experimental  measurements)  may  be  obtained  from  the  University  of  Houston's 
Data  Bank,    compiled  by  Dukler  et  al.    (15). 


Levy'  s  approximation  theory  may  be  found  on  page    54. 


ANALYSIS  (TWO-PHASE  FLOW) 

Before  analyzing  two -phase  pressure -drop  correlations,    it  will  be 
necessary  to  describe  the  various  flow  classifications  which  are  most  gener- 
ally accepted  in  the  literature.     These  classifications  provide  the  means 
whereby  mathematical  expressions  may  be  applied  to  the  energy  and  momen- 
tum balances.     Two-phase  flow  is  classified  by  system,    appearance,    and 
type. 

System  Classifications 

Classification  by  system  involves  either  being  single  component  (a  pure 
liquid  and  its  vapor),    or  two  or  more  components  with  any  one  component 
present  in  either  or  both  of  its  phases.     Systems  may  have  mass  transfer 
between  phases  (i.  e.,    vaporization),    or  otherwise.     They  may  be  isothermal, 
adiabatic,    or  have  intermediate  temperature  behaviors. 

Appearance  Classification 

A  flow  pattern  refers  to  the  visual  appearance  of  a  particular  flow  con- 
dition.    Alves  (2)  illustrated  seven  flow  patterns  for  horizontal  flow  (Fig.  1). 
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Fig.  1.     Flow  patterns  in  horizontal  flow. 


These  flow  patterns  are: 

Bubble  flow.      Discrete  gas  bubbles  move  along  the  upper  surface  of  the 
pipe  at  approximately  the  same  velocity  as  the  liquid.     At  high  liquid  rates, 
bubbles  may  be  dispersed  throughout  the  liquid,  a  pattern  often  referred  to  as 
froth  flow . 

Plug  flow.      The  gas  bubbles  tend  to  coalesce  as  gas  flow  rate  increases, 
to  form  gas  plugs  which  may  fill  a  large  part  of  the  cross -sectional  area  of 
smaller  tubes. 

Stratified  flow.     Complete  stratification  of  gas  and  liquid,  with  the  gas 
occupying  a  constant  fraction  of  the  cross-sectional  area  in  the  upper  portion 
of  pipe,  over  a  smooth  liquid-gas  interface.     It  occurs  at  lower  liquid  rates 
than  bubble  or  plug  flow,  and  more  readily  in  larger  tubes. 

Wavy  flow.    Increasing  gas  rate  produces  waves  of  increasing  amplitude 
at  the  stratified  gas-liquid  interface,  because  of  the  higher  gas  velocity. 

Slug  flow.    Wave  amplitudes  increase  to  seal  the  tube,  and  the  liquid 
wave  is  picked  up  by  the  rapidly  moving  gas  to  form  a  frothy  slug  which  passes 
through  the  pipe  at  a  much  greater  velocity  than  the  average  liquid  velocity. 
Slug  flow  is  also  formed  from  plug  flow  as  the  gas  flow  rate  is  increased  at 
constant  liquid  rate. 

Annular  flow.     Gravitational  forces  become  less  important  than  inter- 
phase forces,  and  the  liquid  is  mainly  carried  as  a  thin  film  along  the  tube 
wall.     The  gas  moves  at  a  high  velocity  in  the  core  of  the  tube  and  carries 
with  it  some  of  the  liquid  as  a  spray.     Film  flow  is  a  name  also  applied  to 
this  pattern. 

Mist  or  spray  flow.     More  and  more  liquid  is  carried  in  the  gaseous  core 
at  the  expense  of  the  annular  film,  until  nearly  all  of  the  liquid  is  entrained  in 
the  gas.     The  pattern  has  also  been  called  dispersed  flow  or  fog  flow. 


10 


Upward  vertical  flow  patterns  are  more  complex  than  horizontal  flow 
patterns  because  as  the  gas  velocity  increases  at  a  constant  rate,  a  dispersed 
type  of  flow  will  be  reached  at  lower  gas  velocities  due  to  the  influence  of  gravity 
on  the  liquid.     Vertical  flow  patterns  tend  toward  radial  symmetry.     This  is 
not  the  case  in  horizontal  flow.     A  classification  of  vertical  flow  patterns 
based  largely  on  air-water  mixtures  has  been  given  by  Nicklin  and  Davidson 
(39)  (see  Fig.   2).       (No  literature  on  downward  flow  was  found.) 
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Fig.  2.     Flow  patterns  in  vertical  flow. 


Bubble  flow.     Gas  is  dispersed  in  the  upward  flowing  liquid  in  the  form 
of  individual  bubbles  of  various  sizes.     As  gas  flow  increases,  the  bubbles  in- 
crease steadily  in  numbers  and  size. 

Slug  (or  plug)  flow.      The  gas  bubbles  coalesce  to  form  larger  bullet 
shaped  slugs  having  a  parabolic  outline  at  the  head.     These  slugs  increase  in 
length  and  diameter,   and  their  upward  velocity  increases  as  the  gas  rate  in- 
creases.    The  slugs  are  separated  by  liquid  plugs  which  contain  gas  bubble 
inclusions.     As  the  gas  slug  moves  along  the  tube,   liquid  flows  down  through 
the  thin  liquid  annulus  surrounding  it  into  I  lie  bubble-filled  liquid  plug  beneath, 


11 


Froth  flow.      When  back  flow  of  the  liquid  around  the  slugs  nearly  stops 
the  slug  becomes  unstable,   and  gas  slugs  seem  to  merge  with  the  liquid  into  a 
patternless  turbulent  mixture  having  the  general  nature  of  a  coarse  emulsion. 
The  elements  of  this  structure  are  in  a  continual  process  of  collapse  and  re- 
formation. 

Annular  (or  climbing  film)  flow.    The  gas  travels  up  the  core  of  the  tube 
at  a  high  velocity  and  the  liquid  forms  an  annular  film  around  the  tube  walls. 
Initially,  this  film  may  be  fairly  thick  and  have  long  waves  on  which  are 
superimposed  a  pattern  of  fine  capillary  waves.     As  gas  flow  rate  increases, 
the  film  becomes  thinner  and  the  amount  of  the  liquid  entrained  as  droplets 
in  the  central  gaseous  core  increases. 

Mist  (spray,  fog,   fully  dispersed)  flow.     At  very  high  gas  rates  the 
amount  of  liquid  entrainment  increases  until  apparently  all  the  liquid  is 
carried  up  the  tube  as  a  mist  or  fog.     Although  a  thin  liquid  film  may  exist  on 
the  wall,   its  presence  is  not  obvious  in  this  region. 

Distinguishing  the  various  possible  flow  patterns  has  usually  been 
approached  experimentally  by  visual  observation  where  the  results  represent 
the  space  relationships  of  the  two  phases.     Photographs,  moving  pictures, 
and  stroboscopic  techniques  have  been  utilized  to  more  clearly  define  transi- 
tions between  the  various  flow  patterns.     So  far  these  transitions  have  not 
been  acceptably  defined. 

Type  Classification 

Because  of  the  inherent  difficulties  present  in  classifying  and  applying 
flow  patterns  to  a  correlation  model,    Martinclli  et  al.    (35)  classified  two- 
phase  flow  on  the  basis  of  turbulence^r  lack  of  it^which  exists  in  each  phase. 
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This  classification  results  in  four  possible  combinations  hereafter  referred 
to  as  flow  types .       They  are: 

1.  Viscous  liquid  flow  —viscous  gas  flow  (v-v). 

2.  Viscous  liquid  flow  —turbulent  gas  flow  (v-t). 

3.  Turbulent  liquid  flow  —  viscous  gas  flow  (t-v). 

4.  Turbulent  liquid  flow  —  turbulent  gas  flow  (t-t). 

Lockhart  and  Martinelli  (33)  used  the  Reynold'  s  number  to  determine  the  tur- 
bulence or  non-turbulence  of  each  phase.     (Re  >  2000     -    turbulent;  Re  <  1000  - 
viscous,  where  1000<Re<2000  indicates  transitional  instability.) 

Other  methods  have  been  proposed  to  predict  flow  types  from  flow  con- 
ditions, fluid  properties  and  pipe  geometries.     These  have  been  presented  by 
Baker  (7),  Griffith  and  Wallis  (21),  Govier    et  al.    (19)  among  others. 
From  these  charts  approximate  determinations  can  be  made  as  to  the 
type  of  flow  which  should  occur  for  a  specific  flow  problem.     This  information 
may  then  be  used  in  one  of  two  ways:      either  to  facilitate  physical  comprehension 
of  the    fluid  motion  in  the  equipment,  or  to  allow  the  selection  of  the  best 
pressure-drop  correlation  within  the  various  flow  conditions.     A  sample  prob- 
lem is  included  in  Appendix  A  which  demonstrates  Baker'  s  method.     Anderson 
and  Russell  (3)  present  an  excellent  review  of  this  subject. 

Recently,  Smissaert  (45)  presented  a  paper  in  which  he  stated,       Flow 
pattern  effects  were  found  to  be  conveniently  described  by  means  of  a  variable 
exponent  of  the  Froude  number.     This  exponent  was  observed  to  be  a  linearly 
decreasing  function  of  the  logarithm  of  the  ratio  of  the  volumetric  flow  rates. 

This  experimental  study  applies  only  to  vapor-liquid,    low-circulation 
rates,  but  it  is  indicative  of  the  type  of  research  and  conclusions  that  are 
being  presented  in  an  effort  to  comprehend  the  mechanisms  and  phenomena 
of  the  two -phase  flow. 
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Special  Terms 

In  conjunction  with  the  various  classifications  of  two-phase  flow  pre- 
viously presented,  it  is  useful  to  define  and  derive  some  terms  pertinent  to 
the  study  of  two-phase  flow  and  two-phase  pressure  drop.     (For  description 
of  nomenclature  see  Appendix  D.  ) 

Flow  pattern.      Flow  classification  based  on  the  visual  appearance  of 
the  flow. 

Flow  type .  Distinguishing  between  the  various  flow  mechanisms  on 
the  basis  of  turbulence  or  non-turbulence,  using  Reynold1  s  number  as  the 
criterion. 

Flow  model.     A  theoretical  combination  of  flow  patterns  and  types 
which  enables  utilization  of  simplified  mathematical  techniques  to  solve  for 
unknown  flow  parameters. 

Mixture  quality.     Ratio  of  the  mass  flow  rate  of  the  gas  to  the  total 
mass  flow  rate  of  the  mixture: 
Vapor  mass  fraction    =    Quality    =    X    =    w    /w  (1) 

O 

where  w    =    w      +    w„  (2) 

g  S. 

Void  fraction.      Ratio  of  the  volume  occupied  by  the  gas  to  the  total 

volume  of  a  channel  section  which  brackets  the  cross  section  of  interest. 

The  void  fraction  is  also  given  by  the  ratio  of  the  cross  section  through 

which  the  gas  is  flowing  to  the  total  cross  section  of  the  channel: 

Void  fraction  =  a   =  v    /v  =  A    /A  (3) 

g'  g'     P 

Slip  ratio.      The  slip  ratio  at  a  particular  cross  section  is  defined  as  the 
ratio  of  the  actual  gas  velocity  and  the  actual  liquid  velocity  at  the  cross-section: 

Slip  ratio  =  u    In  .  (4) 

g       f 

Volume  fraction.        Fraction  of  pipe  occupied  by  a  single  phase  of  a  mixture. 
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Holdup.      Holdup  predicts   the   fraction  of  the  pipe  occupied  by  a   single 
phase  of  the  mixture  at  a  specified  cross   section  in  the  pipe  for  known  quantities 
of  liquid  and  gas  entering. 

Relative  velocity.     Defined  at  a  particular  cross  section  of  the  channel, 

relative  velocity  is  the  difference  between  the  actual  gas  velocity  and  the 

actual  liquid  velocity  at  the  cross  section: 

Relative  velocity  =  u'    =    u     -   u„  (5) 

Superficial  liquid  velocity.    The  liquid  velocity  calculated  on  the  basis 

of  the  liquid  mass  flow  rate,  the  liquid  density,   and  the  total  cross  section  of 

the  conduit: 

Superficial  liquid  velocity  =  u     =  w. /A  p  (6) 

S  x         ] )    Si 

Superficial  velocity  of  the  mixture.     Mixture  velocity  calculated  on 
the  basis  of  the  total  mass  flow  rate,  the  liquid  density,   and  the  total  cross 
section  of  the  channel.     (This  rather  strange  combination  of  terms  finds  its 
origin  in  the  fact  that  w     is  generally  a  much  smaller  quantity  than  w     .     The 
superficial  velocity  of  the  mixture  is,  in  most  cases,   approximated  by  the 
superficial  liquid  velocity.      This  approximation  does  not  hold  for  low  circu- 
lation rates. ) 

Superficial  velocity  of  the  mixture  =  u         =  w/A    n„        (7) 
i  J  sm  '     pKi>        v    ' 

The  superficial  gas  velocity.      The  superficial  gas  velocity  at  a 

particular  cross  section  is  the  gas  velocity  calculated  on  the  basis  of  the  mass 

flow  rate  of  the  gas,  the  total  channel  cross  section,  and  the  density  of  the 

gas  at  the  cross  section: 

Superficial  gas  velocity  =  u        =  w    /p   A  (8) 

f  *  J         sg  g/Kg    p 

Derivation  of  Fundamental  Equations 

Because,   in  various  two-phase  pressure  drop  methods,   the  investiga- 
tion will  encounter  equations  involving  holdup,   slip  and  relative  velocity,   it 
is  important  to  recognize  the  different  forms  by  which  these  terms  are 
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expressed.     The  following  summary  will  be  helpful.     The  gas  and  liquid 
velocities  are  defined  to  be: 


and 


where 


and 


u     =  w    /A  p 

g         g      gg 


U  „  =    W  J  An 


A      =  a  A 
g  P 


(9) 


(10) 


(11) 


Ai-(l-«)Ap. 

The  slip  ratio  is,  therefore,  equal  to 

w      . 

u    /u    =  — &  

g     i      w  a 
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The  ratio  of  the  mass  flow  rates  can  also  be  expressed  as 


w    /w    = 
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w    /w 
SI 
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g'     $_      w     -    w  1     -     (w    /w)  1 
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and  the  slip  ratio  is  then  determined  by 
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(12) 
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(14) 


JS    =  X  1    -   a      H 

Un  1      -     X  n  p 
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(15) 


Subtraction  of  Eq.    (10)  from  Eq.    (9)  and  substitution  of  Eqs.    (11)  and  (12) 
give  the  relative  velocity 
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u     -  u 


g  o      a  A   p 

P    g 


(1     -   a)  ApPi 


(16) 
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u     -  u„  =  w 

g     4 


w                 , 
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Substitution  of  Eq.    (14)  and  simplification  gives 
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u      -  u  „  = 


g        i     ApPf 
Substitution  of  Eq.    (7)  yields 
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(17) 


(18) 


(19) 


16 


Physical  Relationships  Among  Flow  Parameters 
To  acquaint  the  reader  with  a  physical  picture  of  two-phase  flow  inter- 
action, the  following  description  given  by  Smissaert  (45)  outlines  the  funda- 
mental relationships  among  the  various    flow  parameters: 
1 .      Influence  of  Flow  Conditions 

As  a  general  rule,   the  slip  ratio  and  the  relative  velocity  increase 
with  increasing  quality.     However,  the  rate  of  change  of  the  slip  as  a 
function  of  quality  is  larger  for  low  qualities  than  for  high  qualities.     The 
influence  of  quality  on  the  slip  ratio  decreases  significantly  with  increasing 
pressure  and  also  decreases  with  increasing  circulation  rate,  but  to  a  lesser  exten 

'Slip  ratios  decrease  and  relative  velocities  increase  with  increas- 
ing superficial  liquid  velocity.     The  magnitude  of  the  velocity  effect 
appears  to  decrease  with  decreasing  quality  and  increasing  pressure. 
The  rate  of  change  of  the  slip  as  a  function  of  quality  also  decreases 
with  increasing  circulation  rate. 

Pressure  has  a  dual  influence  on  two-phase  flow  phenomena.     The 
increase  in  gas  density,  which  is  a  direct  effect  of  the  increase  in 
pressure,  will  reduce  the  buoyant  force  which  the  liquid  is  exerting  on 
the  gaseous  phase.     As  a  result,  the  slip  ratio  will  decrease  with  increas- 
ing pressure.     A  second  effect  is  the  decrease  of  the  volume  of  gas 
which  is  present  in  the  two-phase  mixture.     This  effect  results  in  a 
change  of  flow  pattern,  which  has  a  substantial  influence  when  a  churn 
or  semi-annular  flow  pattern  is  effective. 

Geometry  includes  such  items  as  the  form  of  the  channel  cross- 
section  (rectangular,   circular',  etc.  )  and  the  magnitude  of  the  equiv- 
alent diameters  smaller  lhan  2  in.     The  effect  of  geometry  is  more 
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significant  the  smaller  the  equivalent  diameter.     Flow  patterns  are 
believed  to  be  responsible  for  this  phenomenon. 
2 .      Influence  of  the  Fluid  Properties 

Although  previous  studies  have  been  successful  in  establishing  an 
acceptable  theory  on  the  effects  of  flow  conditions,  the  same  cannot  be 
said  about  the  understanding  of  the  influence  of  fluid  properties.     Al- 
though these  effects  have  been  studied  in  the  past,  there  still  exist 
numerous  controversies  about  their  relative  importance.     Fluid  prop- 
erties (density  not  included)  are  considered  to  have  little  or  no  influence 
on  slippage.     It  should  be  noted,  however,  that  the  experimental  study 
of  the  influence  of  fluid  properties  is  inherently  difficult.     There  is 
virtually  no  way  to  change  one  property  of  a  fluid  and  keep  the  others 
constant  at  the  same  time.     Moreover,  in  order  to  investigate  a  suffi- 
ciently wide  range,  one  is  compelled  to  employ  different  fluids.    As  a 
result,  the  influence  of  one  property  is  generally  masked  by  that  of  the 
other  properties. 

"it  is  an  easily  accepted  and  experimentally  verified  fact  that  the 
slip  ratios  increase  with  increasing  liquid  density.     From  a  previous 
study  by  Moore  and  Wilde  (38)  it  seems  that  slippage  is  to  some  extent 
dependent  upon  the  surface  tension  of  the  liquid.     The  influence  of 
viscosity  has  repeatedly  been  reported  as  negligible." 
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ANALYSIS  (TWO-PHASE  PRESSURE  DROP) 
Generalized  Analysis 

By  interaction  between  phases,    two-phase  pressure  drops  are  generally 

greater  than  those  experienced  in  single -phase  flows  with  comparable  mass 

flow  rates.     For  any  tube,    two-phase  pressure  drop  consists  of  the  summation 

of  losses  due  to  friction,   acceleration,    and  elevation.     Frictional  losses  (the 

amount  of  energy  that  the  medium  releases  in  overcoming  friction  between  the 

tube  wall  and  the  fluid)  are  always  present  and  occur  for  any  orientation  of 

the  flow  channel  and  for  adiabatic  or  heated  conditions.     These  losses  are  the 

most  difficult  to  analyze.     Acceleration  losses  (change  of  kinetic  energy 

caused  when  the  medium  accelerates  or  changes  momentum)  which  for  horizontal 

pipes  and  no  friction  usually  are  evaluated  simply  as  A    AP     (momentum  of 

P         a 

inlet  streams  minus  the  momentum  of  outlet  streams).     Momentum  losses 
frequently  occur  in  forced-circulation  boiling  when  vaporization  causes  the 
mixture  density  to  vary  along  the  tube.     Losses  in  pressure  due  to  elevation 
are  associated  with  the  weight  of  the  medium  whenever  it  changes  its  vertical 
elevation  and  are  present  both  for  vertical  and  inclined  flow  systems. 


Allowances  for  hydrostatic  head  also  are  usually  made  separately  by  AP. 

r»  2 
=  -=*-      \       dz/ v  where 
g, 


r 2 

dz/ 
J  1 


1/v  =  PgRg+Pi  <*    "Rg>  <20> 

To  predict  acceleration  losses  and  elevation  losses,    knowledge  of  the 
mean  two-phase  mixture  density  is  necessary.     In  vapor-liquid  flow,    the  mean 
velocities  of  the  two  phases  (dependent  on  the  cross -sectional  area  of  each) 
generally  are  not  equal.     Because  of  this  "slip,"  the  true  frnrtion  of  the  pipe 
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cross -section  occupied  by  either  phase  differs  from  the  ratio  of  gas  and 
liquid  volumes  entering  the  tube.     Asa  result,    the  mean  mixture  density 
cannot  be  calculated  on  the  basis  of  quality  (vapor  mass  fraction)  alone,   but 
requires  a  knowledge  of  the  void  fraction  (fraction  of  the  channel  cross  section 
occupied  by  vapor).     Void  fraction,    then,    is  an  important  parameter  used  to 
predict  the  hydraulics  of  many  two -phase  flow  systems,    since  without  it,   the 
accelerative  and  hydrostatic  contributions  to  the  pressure  drop  cannot  be 
evaluated. 

Holdup  correlations,    to  predict  the  volume  fraction  of  liquid  and  gas  at 
a  specified  cross  section  in  the  pipe  for  known  quantities  of  liquid  and  gas 
entering  the  pipe,    are  essential  for  most  frictional  pressure-drop  calculations. 
Holdup  differs  from  the  volume  fraction  determined  from  the  densities  and 
mass  flow  rate  because  of  slip  and  density  changes  caused  by  the  pressure 
gradient  in  the  pipe. 

Early  methods  to  determine  holdup  suggested  by  Lockhart  and 
Martinelli  (33),    Levy  (32),    Martinelli  and  Nelson  (36),    and  Hughmark  (26)  con- 
sider the  extreme  flow  models  of  homogeneous  and  annular  flow.     Bankoff  (8) 
in  1960  suggested  a  holdup  correlation  which  did  not  assume  these  extreme 
conditions  of  flow.     He  considered  a  model  in  which  the  mixture  flows  as  a 
suspension  of  bubbles  in  the  liquid,    and  where  radial  gradients  exist  in  the 
concentration  of  bubbles. 

By  assuming  a  power  law  distribution  for  the  velocity  and  void  fraction 

u/u       =S1/m  (21) 

1     m 

<*/om  =  Sl/n  (22) 
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Bankoff  (8)  derives  the  following  equation  for  holdup: 


I      1      ui 
x  p 


i 


K 


R 


(23) 


g  / 


The  flow  parameter  K"  is  a  function  of  m  and  n,  and  was  calculated  from  the 
empirical  equation 


K"    =  0.71  +  0.0001  P 


(24) 


K"    applies   only  for    steam-water     (atmospheric  -  2000  lb/in.    )  systems.     This 
equation  does  not  apply  to  air-liquid, two-phase  flow. 

Hughmark  (25)  correlated  a  more  general  method  for  determining  K"  in 
Eq.  23.   His  flow  parameter  (designated  by  K)  is  obtained  from  the  dimensionless 
chart  shown  in  Fig.  3  which  relates  the  following  three  variables: 
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Fig.  3.    Hughmark' s  correlation  for  flow  parameter  K  from  Z. 
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1.     Reynold's  Number 

ReTP  ■  (iVVTyT)  (25) 


where  (R,,M,>  +  Rju    )  is  a  two-phase  viscosity  term. 

2.  Froude  Number 

2 
tt  UNS  ,„   t 

Frm  =  Jd"  (26) 

(the  velocity  is  defined  for  the  two  phases  traveling  at  the  same  velocity, 
i.  e.,   no  slip  between  phases). 

3.  Liquid  volume  fraction 


w   v 

0      0 

X    =  ^~ ,  (27) 

w„    v.   +  W    V  x       ' 

a   i       g  g 


by  the  equation 


z  =  1? Hi .  (28) 

oo1/4 


To  determine  the  volume  fraction  (holdup),    an  example  problem  is 

12 

included  as  a  part  of  Dukler's  method  to  predict  pressure  drop. 

The  experimental  techniques  currently  employed  to  measure  pressure 
drop  in    two-phase  flow  involve  the  use  of  piezometer  rings  for  vertical  flow, 
and  single -pressure  taps  for  horizontal  flow.     Normally,   phase  separators 
adjacent  to  the   pressure  tap  ensure  single -phase  communication  with  the  pres- 
sure measuring  device.     Essentially,   the  techniques  involved  are  strongly 
dependent  on  the  flow  pattern  under  investigation.     Special  problems  occur 


1 2 

See  page  33  for  Dukler's   method  and  example  problem. 
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when  flow  patterns  of  the   plug  and  slug  type  are  considered.     The  mean  pres- 
sure at  any  cross  section  fluctuates  considerably,    and  the  determination  and 
significance  of  a  mean  pressure  with  time  poses  difficult  problems.     For 
steady  horizontal  flow  (i.  e.    stratified  flow)  another  difficulty  arises  due  to 
the  variation  in  pressure  around  any  vertical  diagnostic  tube.     Because  it 
is  only  in  annular  and  dispersed  flow  that  these  difficulties  are  reduced  to 
a  minimum,    annular  flow  has  received  the  most  experimental  and  theoretical 
study. 

A  comprehensive  listing  of  recent  experimental  steam-water  investi- 

1 3 
gations  has  been  given  by  Tong  (48).         A  similar  table  has  been  presented 

by  Hughmark  (25),    including  different  two-phase  components. 

Energy  and  Momentum  Equations 

Energy  and  momentum  equations  are  fundamental  to  the  analysis, 
derivation,    and  discussion  of  any  two-phase  pressure  drop  calculation  method. 
In  1962,    Vohr  (49)  published  a  paper  in  which  the  proper  form  of  the  energy  equation 
(no  mass  transfer)  is  discussed.    White  and  Lamb  (50)  in  the  same  journal  reviewed 
the  use  of  momentum  and  energy  equations  in  two-phase  flow.     Scott  (43) 
presented  an  excellent  review  in  which  the  limits  and  common  conditions  en- 
countered experimentally  are  briefly  summarized.     This  paper  includes  a  re- 
view of  a  paper  by  Nicklin  and  Davidson  (39)  in  which  the  irreversible  terms 
of  the  equation  are  discussed.     These  articles  furnish  ample  background  for 
applying  the  equations  of  energy  and  momentum  to  the  pressure-drop  cor- 
relations presented  in  the  literature. 

The  energy  equation,    mechanical  energy  equation,    and  momentum  equation 
essential  for  the  analysis,    derivation  and  discussion  of  two-phase  pressure 
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See  Appendix  B. 
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drop  correlation  methods  are  derived  in  a  manner  dependent  upon  whether 
the  flow  pattern  being  considered  is  continuous  (annular)  or  discontinuous 
(plug  flow).     Continuous  flow  patterns  are  analyzed  microscopically. r  Dis- 
continuous  flow  patterns  are  analyzed  with  sufficient  accuracy  by  macroscopic 
balances.        In  either  case,    assumptions  are  usually  made  to  simplify  the 
mathematics  as  much  as  possible.     Since  most  two-phase  flow  investigations 
are  carried  out  with  one  dimension  in  the  steady  state  and  with  constant  flow 
rates,   these  assumptions  don't  significantly  affect  the  accuracy  of  the  resulting 
equations.     For  two-phase  flow  systems,   with  or  without  mass  transfer  be- 
tween phases,      and  assuming  no  shaft  work,    Eqs.    2  9-32  result. 

Macroscopic  Balance 
Macroscopically,   the  usual  approach  assumes  pressure  constant  for  any 
given  cross  section  of  the  pipe.     Momentum  and  energy  equations  are  then 
written  separately  for  each  phase  with  the  constraint  that  the  static  pressure 
drop  is  identical  for  both  phases  over  the  same  increment  of  flow  length,   then 
added  to  give  overall  expressions.     Unfortunately,    the  resulting  overall  two- 
phase  balances  do  not  have  simple  relationships  to  each  other  as  exist  in  single - 
phase  flow.     It  becomes  necessary  to  weight  the  momentum  terms  by  their 
respective  velocities  in  order  to  form  correct  energy  quantities.     By  adding 
these  weighted  equations  which  are  applicable  to  each     phase,    the  results  are: 

a.     Momentum  Equation 

f  \     d(wiV     d(wgV  .  e 

-  AdP  =  -    A.dP„,    +AdPu,       = Z-JL  + £_JL_   +    p   A .  +  p   A  )  -*-   dz 

V  s.     w        g     w  )        g  g  Vfjc  i     ' g  g7  g 


+  <Tw   c'w  +Yw   cw    'dz  (29) 

g        g        i       i 


See  footnote,    p.    24. 
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b.     Mechanical  Energy  Equation 


-    2' 


d(Vg)  +  d(Vui  ) 


g  g 


2g, 


+  (A„u      +  A    u    )  dP  +  (p.A.u      +  p   A    u    ) 


-&-     dz 


-   -dEg  -  dE^   =   -dEgi  (30) 


c.     Energy  Equation  (no  shaft  work,    adiabatic) 


_d_ 
dz 


=  0    (31) 


After  performing  the  differentiation  (noting  that  z     =  z       dp.   =  0, 

and  dw    =  -dwfl),    there  results: 
g  & 


[  A    u    P  ]  — &  +  w    dU     +  w„  dU„  +     &   A — £-     + 
I    g   g     )   P0  g       g         i       i 


/u2         u2  > 

+  (w     +w„)f  dz+f- -")  Pdw     +(u     -U„ldw     +    o^-  -  t^- 

g         *     gc  \Pg      Pjt)  g       V    g         i^       g      \^2gc       2gcj 


dw     =  0 
g 


(32) 


Microscopic"  Balance 
Microscopically,    continuous  flow  may  be  analyzed  by    differential  equations 
simplified  according  to  the  various  simplifying  assumptions  applicable  to 
the  specific  problem.      The  following  example  follows  the  unpublished  notes 

of  Mr.    G.    W.    Patraw: 

For  a  system  having  constant  area  flow  with  phase  change,  friction, 
heat  addition,  and  considering  a  pipe  element  dL,  the  inlet  parameters  to 
a  fluid  element  of  unit  width  and  length  dL  are  given  in  Table  1. 


Microscopic  balance  refers  to  stratified,    annular,    and  spray  horizontal 
flow  patterns  (Fig.    1)  as  well  as  annular  and  mist  vertical  flow  patterns  (Fig.    2). 
Macroscopic  balances  refer  to  the  remainder  of  the  flow  patterns. 
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Table  1.     Inlet  parameters  to  a  fluid  element. 
Property Vapor  Liquid  Total 

wihi>  Wghg  +  Wi!hi 


vi 


Thermal  equilibrium  between  the  two  phases  is  assumed. 
1.     Energy  Balance  Across  the  Element 
a.     Differential  change  in  enthalpy 

dH  =  d  (w^hp)  +  d  (w    h  ) 

dH  =  h^  dw^  +  Wje  dh^  +  hg  dwg  +  wg  dhg  (33) 


Velocity 

u 
g 

Mass  Flow 

wg 

Enthalpy  (static) 

w    h 
g    g 

Pressure  (static) 

P 

Specific  Volume 

V 

b.  Differential  change  in  kinetic  energy 

dKE  =  5?7  d K  »ia  +  Wg \) 

dKE  =  sir  W  dwi +  2  wi  ui  du* +  ug2  dwg +  2  wg d  ug)         (34) 

c.  Heat  added  to  system 

dQ'  =  qdL  (35) 

There  are  no  potential  energy  changes  or  shaft  work  so 

dH  +  dKE  =  dQ'  (36) 

(  w   ^ 

Introducing  continuity  and  quality  (  X  =  — *»J  (1) 

w     =  Xw  dw     =  w  dX  (la) 

g  g 

w     =  (1   -  X)w  dw^  =  -w  dX  (lb) 

and  letting  h  „  =h     -  h„,    the  heat  of  vaporization,   the  result  is:  (37) 

gJ?        g         * 
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w 


h    +-^-r    fu  2  -  un2)    dX 
_g^    2gcJ     V    g  i  JJ 


+  w(l    -  X)  dh„  +  wX  dh 
i  g 


+    W 


(1   -  X)  u„  du     +  X  u     du 


q  dL=  0    (38) 


gc.T     [_  '        JL        I  g       g 

For  most  flows  of  interest,    the  kinetic  energy  terms  are  very  small  com- 
pared to  the   thermal  energy  terms,    so  the  energy  equation  becomes 

w  h    „  dX  +  w(l   -X)  dh     +  wX  dh     -  q  dL=  0  (39) 

gi  i  g 

The  change  in  enthalpy  terms  (dh)  are  also  small,    so  we  have  approximately 
on  integration 

Q'    =  wh    „  (X„.      .   -  X.    ...    ,)  (40) 

gi       final         initial 

which  is  the  usual  approximation  made  in  a  heat  flow  balance.     For  this 

particular  analysis,    the  energy  equation  need  no  longer  be  considered,    other 

than  remembering  it  shows  that  the  quality  is  linear  with  distance  for  constant  heat 

input  (  -Tf-  =  constant  ]  . 

2.     Momentum  Balance  Across  the  Element 

The  change  in  pressure  is  equal  to  the  change  in  pressure  due  to  friction 

plus  the  change  in  pressure  due  to  the  change  in  momentum  of  the  stream.     (dPf 

is  always  negative.) 

dPTp  =  dPf+dPa  (41) 

The  dPf  term  will  be  discussed  later.     The  momentum  term  is: 

-dP       = — ~—    (1   -  X)  w  dufl  +  Xw  du     +  (u     -  uj  wdX  (42) 

a        gcAp    (_  JL  g  g        i  J 

The  liquid  and  gas  phases  may  flow  at  different  velocities.     From  con- 
tinuity 

(I   -  X)  v.  w  X  v  w 

ut- — ^ —        V"*f"  <43,44) 
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Substituting  Eq.    (11)  into  Eqs.    43,    44  we  obtain 


u 


(1    -  X)  Vje  w 


1       (1   -  a)  A 


P 


du„   = 


(1   -  X)  v^  w 


1       (1   -a)  A. 


u     = 
g 


X  v    w 

g__ 

or  A 


dv. 


dX 


+ 


eta 


1   -  X        1   - 


a 


(45,46) 


(47) 


du     =-S- 


v    wX    /dv 


g      a  Ar 


+ 


g 


dX         do- 
X  a 


(48) 


-dP     = 
a 


w 


g  A  * 

6C     p  L_ 


(1   -  X)"  vi 
(1   -a) 


dv. 


dX 


+ 


eta 


X     v 


+ 


A 


dv 


£  +  ^X   _  d^ 


1   -  X        I  -a 

Xv  (1  -X)  v, 


a  X 


+ 


A.  _ 


dX 


( 1   ~  a) 

Since  the  compressibility  of  the  liquid  is  generally  small,    the  pressure - 
momentum  equation  becomes 


(49) 


dP  =  dP 


w 


f      gA2 
sc    p     L 


(1-x)  V 


eta 


dX 


(1   -a)      ^(1   -a)  "  (1   -  X) 

rXvg_(l-X)vi 


\    x  v     /  dv      ,v    .y 

g_   I  g.  +  dX      dor  1 


v 


g 


X       a 


letting 


a 


(1   -a) 


dX 


dv         dv  JT1 

_£  =  __£       dp 


dL 


dP      dL 


(50) 


(51) 


TP 


and  converting  to  derivative  form 


dP 


dP, 


w 


dLTp       dL 


g  A 

&C    p    L 


(i  -  XT  vje 


(1   -or)  1 


eta 


dX 


a 


dL        1   -  X    dL 


X     v      /    ,    dv  ,_ 
g_     _1   _^_dP_  +  1 

a       \g    dP    dLTP 


,      ,v        ,      ,     A    (Xv         (1   -  X)vn\ 
JL    <*X    _  I    eta    I  +|        g  _  i  I 

X    dL     "  o    dL  |    U  (I   -  a)  j 


dX 
dL 


(52) 
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thus 


V        g 


2      v2dv    \     ,_,  dP-        2 

w        X         g\   dP  _J      w 

"  dL  " 


A2  a      dPi  dL 
c    p  '         TP 


gA2 


dLl  (I-.)2 


X' 


a 


2    Vg 


.„  '  Xvff        (1   -  X   v 
dL  I    a 


(53) 


(1   -a) 

This  is  a  form  of  the  momentum  equation  which  may  be  applied  to  the  Homogeneous 
model  presented  on  page  40  of  this  report. 
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CORRELATION  METHODS 

General 

For  every  problem  which  has  resisted  definitive  solution  over  a  long 
period  of  time,    the  methods  of  approach  are  usually  sophisticated  by  steps  in 
which:     (1)  imperfect  correlations  are  presented  based  on  limited  data;    (2) 
these  correlations  are  improved  as  the  problem  receives  further  study  and  more 
sophisticated  tools  are  used  for  its  analysis;  and  (3)  ultimately,   physical  limits 
upon  the  problem  render  various  solutions  more  satisfactory  than  others,   de- 
pending upon  the  magnitude  of  their  complexity  and  application. 

In  the  interim,    five  classical  patterns  may  be  separately  identified: 

1.  Empirical  correlations; 

2.  Correlations  utilizing  dimensional  analysis; 

3.  Correlations  utilizing  similarity  analysis  and  model  theory; 

4.  Mathematical  analysis  of  a  simplified  physical  model  and  development 
of  equations  relating  the  variables; 

5.  Approximation  theory(solutions  to  the  energy,    momentum,    and  con- 
servation equations  with  empirical  expressions  used  for  the  turbulent  transport 
terms,    approximations  to  the  boundary  conditions,    and  assumptions  involving 
the  relative  magnitude  of  various  terms  in  the  equations.     Usually  the  resulting 
relationship     among   variables  is  obtained  by  numerical  solutions  for  the  complex 
equations). 

So  far,    two -phase  pressure  drop  correlations  have  been  published  which 
involve  four  of  the  above  five  types.     The  fifth  appears  unlikely.     All  five 
approaches  are  discussed.     Examples  have  been  shown  which  demonstrate  the 
most  important  methods. 
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Empirical  Approach 

Many  empirical  correlations  have  appeared.     Most  of  these  can  be  used 
with  poor  reliability  beyond  the  range  of  the  data  from  which  they  were  con- 
structed,   as  shown  in  a  paper  by  Duklor  et  al.    (14).     Empirical  correlations 
for  horizontal  flow  appear  in  references  by  Chenoweth  and  Martin  (11), 
Hoogendorn     (24),   Isbin  et  al.    (27)  and  Sobocinski  and  Huntington  (46)       Correla- 
tions for  vertical  conduits  are  given  in  references  by  Calvert  (10),   Govier  and 
Omer  (18)  and  Hughmark  (25).     Excellent  summaries  may  be  found  in  Tong  (48) 
and  Scott  (43). 

Chenoweth  and  Martin  Correlation 

This  correlation  applies  for  turbulent -turbulent  horizontal  flow  in  pipes 
and  was  intended  to  improve  upon  previous  methods  especially  suited  for 
predicting  high  pressures  in  large-diameter  pipes.     Entirely  empirical,    it 
yields  AP^p/AP,     as  a  function  of  liquid-volume   fraction  of  the  feed  with  a 
quantity  ip  f    p0/ip.p     as  a  parameter.     For  this  correlation  AP   '  is  evaluated  as 

©  **•         **-     o 

the  pressure  drop  based  on  the  total  mass  flow,    using  the  liquid  phase  properties. 
The  parameter  ip     PnW  0  P     is  defined  as 

g         X  X         g 


*JL  Pg         AP*    "LC    z/D  +  SK'  lpg 


(54) 


x  L* 

where  the  f    are  evaluated  at  the  total  mass  flow  from  a  Moody  chart,      and  the 
term  ZK' is  an  allowance  for  valves  and  fittings. 

The  originators  claimed  an  agreement  with  the  experimental  pressure  drops 
of  ±35%  for  this  correlation.     Additional  evaluations  were  made  by  Collier  and 
Hewitt  (12)  who  applied  the  method  to  nir -water,  vertical-annular  flow.     Aziz 
and  Govier  (6)  found  poor  agreemrnt  for  low  liquid  volume  fractions,    and 
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Isbin  et  al.    (27)  investigated  steam -water  flows  at  pressures  from  400-1000 
psia,    and  concluded  errors  over  100%  occurred  over  much  of  the  liquid- 
volume   fraction  range  (predictions  being  low  at  low  volume -fractions,    and  high 
at  high  volume  fractions).    The  net  result  of  these  later  investigations  indicated 
that  this  correlation  gives  results  somewhat  worse  than  the  Lockhart- 
Martinelli  methods  which  are  presented  on  page  43. 

Hughmark  Correlation 

Hughmark  (25)  proposed  a  similar  correlation  which  plots  the  ratio  of 
the  friction  factors  for  two-phase  flow  and  the  flow  of  a  fictitious  single -phase 
fluid  against  the  liquid  volume  fraction  entering  but  uses  a  slightly  different  param- 
eter,   that  is,   (Re  Tp_m)     (P  /P J  (0.085/D)  as  a  third  quantity.     This  choice 
of  parameters  amounts  to  using  a  fluid  viscosity  weighted  for  the  quantities 
of  the  two  fluids  rather  than  a  ratio  of  viscosities  as  in  Chenoweth  and  Martin's 
work,    and  introducing  a  diameter,   inasmuch  as  Hughmark' s  work  was  based 
on  1-in.   pipe.     The  correlation  is  claimed  to  be  generally  applicable  for 
horizontal  flow,   but  in  comparison  with  Chenoweth  and  Martin's  work,   the 
improvement  shown  is  slight. 

Dimensional  Method 

The  general  use  of  dimensionless  groups  for  empirical  correlations  is 
widespread;  however,    two -phase  pressure  drop  correlations  based  primarily 
on    dimensional  analysis  have  not  yet  appeared  in  the  literature.     This  is  due 
to  the  number  of  variables  which  exist  in  a  two-phase  flow.     Among  the  in- 
dependent variables  involved  arc: 
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1.     diameter, 
2,3.     mass  flow  rate  of  each  phase, 
4,5.     density  of  each  phase, 
6,7.     viscosities  of  each  phase, 

8.  interfacial  tension.    ;md 

9.  acceleration  of  gravity. 

For  the  above  variables,    6  dimensionless  ratios  may  be  derived.    Neglecting 
gas  viscosity  and  interfacial  tension,    4  dimensionless  ratios  still  exist. 

When  heat  transfer  is  present  the  thermal  properties  of  each  phase, 
such  as  thermal  conductivity  and  specific  heat,    as  well  as  the  temperature 
profiles  of  the  system,  must  be  considered.    When  mass  transfer  is  taking  place, 
additional  factors  will  include  the  diffusion  coefficient  in  each  phase  for  each 
component  being  transferred,    the  concentration  of  each  component  in  each 
phase,    and  the  equilibrium  relationships  between  phases  for  each  component. 
This  very  large  number  of  factors  explains  why  a  universal  or  general  cor- 
relation for  any  transport  process  would  be  very  complex,    and  why  dimensional 
analysis  is  at  a  disadvantage  in  comparison  with  other  experimental  correlation 
methods. 

Similarity  Analysis  (Dukler  Method) 

Prior  to  a  paper  presented  by  Dukler  et  al.    (15)  in  1964,    no  solution  in 
which  the  principles  of  similarity  in  a  formal  manner  were  used  appeared  in 
the  literature.     An  excellent  review  of  Dukler' s  method  has  been  presented  by 
Anderson  and  Russell  (3). 
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Essentially,    Dukler  relates  the   parameters  involving  two-phase  flow 
by  the  Euler  and  Reynolds  numbers  familiarly  used  in  single -phase  flow. 
Similarity  theory  states  that  two  single-phase  How  systems  are  dynamically 
similar  if  the  Reynolds  numbers  and  the  Euler  numbers  are  equal  respectively. 
In  the  case  of  single -phase  flow,    the  Euler  number  is  twice  the  Fanning  friction 
factor.     This  condition  does  not  provide  enough  information  with  which  to  find 
the  Reynolds  number  or  two -phase  friction  factor.     However,    once  relationships 
are  found  from  experimental  data  for  one  system  (the  model),    the  condition  of 
dynamic  similarity  requires  that  these  same  relationships  apply  to  all  similar 
systems.     In  Dukler1  s  original  paper,    similarity  relationships  for  single-  and 
two-phase  flow  are  developed  in  parallel,  clearly  demonstrating  this  analogy. 

Once  similarity  analysis  has  been  used,    correctly  defining  the  parameters, 
experimental  data  may  then  be  utilized  to  develop  further  relationships  between 
these  parameters.     Thus,    data  from  the  model  may  be  used  to  expand  data 
taken  at  other  conditions  and  in  other  systems.     Dukler  distinguishes  between 
the  various  flow  patterns  by  selecting  constants  and  eliminating  particular  terms 
in  his  original  equations  which  hold  for  the  special  cases  —  no  slip  (homogeneous 
flow),    and  special  slip  conditions.     He  then  statistically  analyzes  the  deviation 
between  theory  and  experimental  values  and  compares  these  results  to  the  method 
of  Lockhart  and  Martinelli  (33).     Dependent  upon  improving  methods  to  predict 
holdup,    Dukler1  s  method  gives  better  agreement  with  experimental  data  than  do 
earlier  correlations  by  Martinelli  et  al.    (35). 

After  the  volume  fraction  has  been  determined,    the  system  of  Dukler  et  al. 
(15)  may  be  used  to  calculate  pressure  drop.     The  following  steps  appear: 
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1.  A  two-phase  Reynolds  number  based  on  flow  conditions,    as  well  as 
the  conditions  in  the  pipe,    must  be  calculated. 

2.  A  friction  factor  is  then  calculated  from  the  two-phase  Reynolds 
number  in  the  Koo  equation  (31). 

3.  The  frictional  pressure  drop  is  calculated  from  an  equation  which 
relates  the  friction  factor,  total  mass  flow  rate,  and  physical  properties  to 
the  frictional  pressure  drop. 

4.  The  effect  of  elevation  is  calculated  by  using  an  average  mixture 
density. 

5.  An  iterative  procedure  is  used  to  evaluate  the  acceleration  correction 
factor.     For  the  case  in  which  the  gas  and  liquid  volume  fractions  and  the  gas 
density  do  not  change  much  along  the  length  of  the  pipe,    the  acceleration 
correction  term  is  negligible.     This  happens  when  the  pressure  drop  is 
small. 

6.  The  individual  contributions  of  friction,    momentum  change,    and  ele- 
vation are  then  added  to  give  the  total  pressure  drop. 

An  example  problem,    presented  by  Anderson  and  Russell  (3),    demonstrates 
Dukler's  method: 

"Given:     a  350-ft  section  of  horizontal  1 -in.    smooth  pipe.     The  flow 
rates  are  as  follows:  1000  lb/hr  of  water  and  15  lb/hr  of  air.     The  discharge 
end  of  the  pipe  is  at  a  pressure  of  1  atm  and  the  system  is  isothermal  at  a 
temperature  of  68°F. 

"Find:     The  predicted  pressure  drop  for  the  test  section  by  using  the 
general  correlation  of  Dukler  et  al.    (14). 
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I.     With  the  method  of  Hughmark  (26)  used  to  calculate  the  volume  fraction 
(holdup)  of  gas  in  the  pipeline: 

'  To  obtain  average  physical  properties  and  to  obtain  an  initial  correction 
for  the  acceleration  losses,    it  is  necessary  to  assume  a  total  pressure  drop 
and  calculate  an  entrance  pressure  and  an  average  pressure.     A  pressure   drop 
of  4.6  psi  will  be  assumed,    giving  an  entrance  pressure  of  19.3  psia  and  an 
average  pressure  of  17  psia.     The  calculations  will  be  done  with  these 
assumptions  and  they  must  be  validated  at  the  end.     The  physical  properties 
at  the  assumed  average  value  for  the  pressure  of  17  psia  are  as  follows: 

p     =  62.4  lb/ ft3,    p     =  0.0870  lb/ ft3,  y     =  1  cP,    and  /u      =  0.018  cP. 

"Assume  the  volume  fractions  in  the  pipe  to  be  R     =  0.75  and  R.   =  0.25. 
Calculation  of  these  fractions  involves  an  iterative  procedure  because  of  the 
dependence  of  Hughmark' s  parameter  Z  on  the  pipeline  volume  fractions.     A 
reasonable  value  for  the  first  assumption,   for  the  liquid-volume  fraction  of  the 
fluid  flowing  in  the  pipe,   is  the  liquid -volume  fraction  of  the  mixture  entering 
the  pipe.     The  calculation  for  the  final  iteration  is  as  follows:    Calculation  of 
the  total  mass  velocity  G: 

G  =  G„  +  G„  (55) 

G  =  1000/(0. 25tt  X  0.08332)  +  15/(0. 25?r  X  0.08332) 
G  =  186,000  lb/hr-ft2 


"Calculation  of  a  two -phase  Reynolds  number: 

ReTp  =  DG/dty^  +  RgMg) 

(0.0 8 3 3)  (186,000) 


(25) 


Re 


TP       (0.25)     (1)  (2.42)  +  (0.75)  (0.018)  (2.42) 


Re         =  35,000 
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"Calculation  of  uNC,>  the  velocity  defined  for  no  slip  between  the  phases: 

UNS  =  <S+Qg>/Ap  (7) 

Q^   =  1000/(62.4  X  3600)  (9) 

Qi    =  0.00445  ft3/ sec 

Q      =  15/(0.087  X  3600)  (10) 

Q      =  0.0479  ft3/  sec 
^g  ' 

uns  b<Q*+Vap 

uN1<3    =  (0.00445  +  0.0479)/(0.25tt  X  0.08332) 

UNS    =  9"60  ft/sec 

"Calculation  of  Fr     ,    the  mixture  Froude  number  based  on  velocity,    which  is  in 

m  J 


turn  calculated  assuming  no  slip  between  the  phases: 

2 

"m      UNS 


Fr       =  u*   /gD  (26) 


Frm  =  9.602/(32.2  X  0.0833) 
Frm   -  34.3 
Calculation  ofX  ,   the  liquid  volume  fraction  of  the  fluid  flowing  in  the  pipe 

evaluated  at  the  average  pressure: 

A  =Qi/(Qg  +  Qi)  (27) 

X  =  0.00445/(0.00445  +  0.0479) 
A  =  0.085 
"Calculation  of  Z: 

Z-(ReTp)1/6(Frm)1/^(Xi)1/4  (28) 

Z  =  (35,000)1/6(34.3)1/8/(0.085)1/4 

Z  =  16.4 
"From  Fig.  3  or  Table  1,    K  is  found  to  be  0.814. 
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"Calculation  of  R    : 

1/X  =  (l   -^j   (1   -K/Rg)  (23) 

X  =  15/1015  =  0.0148 
1/0.0148  =  (1   -  62.4/0.087)  (1    -  0.814/R    ) 


Rg 

=  0.745 

"Calculation 

cfR, 

i  • 

R*  = 

1  -  R 
g 

ni  = 

1   -  0.7 

Rr- 

0.255 

(56) 


From  the  values  obtained  for  R     and  R_,   it  is  seen  that  the  values 

g  a 

assumed  for  these  parameters  in  the  calculation  of  the  two-phase  Reynolds 
number  were  valid. 

"II.     With  the  method  of  Dukler  et  al.    (15)  used  to  calculate  the  frictional 

pressure  drop: 

it 
Calculation  of  P-^c'- 

PNS  =  PjgX  +  Pg(l   -X)  (57) 

PNS  =  (62.4)  (0.085)  +  (0.087)  (0.915) 
PNS  =  5.38  lb/ ft3 
"Calculation  of  /uNC,: 

^NS  =/LiiX  +/ug  (1   "  X)  (58) 

MNS  =  (1)  (0.085)  +  (0.018)  (0.915) 

/uNS  =  0.102  cP 


"Calculation  of  0: 


P  -  (p^/pns)  <x2/iy  +  (pg/pNS>  (1  "X)2/Rg  (59) 

/3  =  (62.4/5.38)  (0.0852/0.255)  +  (0.087/5.38)  (0.9152/ 0.745) 
ft  =  0.360 
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Calculation  of  Re     ,    a  mixture  Reynolds  number: 

m  J 

Rem  -  (DG/^Ng)^ 

Re       =[  (0.0833  X  186,000)/(0. 102  X  2.42)]    [0.360] 

Re       =  32,600 
m 

Calculation  of  friction  factor  f     using  the  Koo  equation  (31): 


(25) 


f     =  0.00140  +  0.125/(Re     ) 
o  '        m 


0.32 


(60) 


f     =  0.00140  +  0.125  /(32.600) 

f     =  0.0059 
o 


0.32 


Calculation  of  F(X)  =  f0/f0  ,    an  empirical  ratio  given  by  Dukler  (15): 

F(X)  =  1.0+7/(1.281  -0.478Y  +  0.444Y2  -  0.044T3  +  0.0084374)     (61) 


where  7  =  -ln(X)  =  -ln(0.085)  =  2.47 

F(X)  =  2.45 
Calculation  of  (AP/AL),.,    the  frictional  pressure  drop: 

(AP/AL)f  =  (2G2fo/gcDpNS]F(X)/3 

/a-o/at  \        11  x  186.0002  X  0.0059  X  2.45  X  0.360) 
(AP/AL).  =  i ' 5 

(32.2  X  3600     X  0.0833  X  5.38) 

(AP/AL)f  =1.91  lb/ ft2  per  ft  of  pipe 

(AP/AL)f  =  0.0133  lb/ in.2  per  ft  of  pipe 
Calculation  of  acceleration  correction  term: 

~~\  r 
(AP/AL). 


(62) 


(63) 


l/gcApAL       wg     Ajl/pgRg]  +    w/-/Pje  Ad/R^)  (64) 


Since  the  assumed  pressure  drop  is  small,    the  changes  R.  and  R     will  be  small, 
and  R     and  R     can  be  taken  as  average  values  of  17  psia.     This  statement  will 
only  be  true  if  the  value  assumed  for  the  pressure  drop  is  the  correct  one.     With 
average  values  for  R.  and  R    are  used,    Eq.    (64)  becomes  the  following: 


(A  P/AL)a  = 
(AP/AL)     = 


1/g  A  'AL 
.       c    p 


i   W  2/R        A(l/p    ) 
g  '     g  ' Mg 


1 
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(65) 


:2//oo    ov     /o,nm2    /n    nr_\2    ,„    nnnnA, 


:  (0. 0833)^(350)  (0.745) 

1         J 
1/0.751   -  1/0.990 

(A  P/AL)      =  0.0000212  lb/ f t^  per  ft  of  pipe 
This  term  is  negligible  compared  to  the  frictional  pressure  drop. 
Calculation  of  (AP/AL): 

(AP/AL)  =  (AP/AL)a  +  (AP/AL)f  (66) 

(AP/AL)  =  0.0133  psi/ft  of  pipe 
Calculation  of  the  actual  pressure  drop: 

(AP)  =  (AP/AL)  AL  (67) 

(AP)  =  (0.0133)  (350) 

(AP)  =  4.65  psi 
The  original  assumption  of  4.6  psi  for  the  pressure  drop  is  valid.     If  it  were 
not,    a  new  value  would  be  assumed  and  the  iterative  procedure  would  be  repeated 
(including  a  new  holdup  calculation)  until  a  satisfactory  agreement  were  obtained. 

Mathematical  Analysis  Involving  Simplified  Physical  Models 

Homogeneous  Model 

Correlations  which  have  received  widest  acclaim  and  popularity  have  been 
those  based  on  a  simplified  physical  model.     The  simplest  mathematical  model 
considers  the  extreme  flow  condition  whereby  the  liquid  and  vapor  phases  are 
completely  and  finely  dispersed.     This  condition  is  approximated  in  the  fully 
dispersed  or  "fog"  flow  pattern.     Due  to  its  simplicity,    this  method  has  been 
applied  to  steam-water  flow  often  without  regard  to  flow  pattern.     The  usual 
assumptions  made  are:     (a)     equal  vapor  and  liquid  velocities;    (b)  attainment 
of  thermodynamic  or  physical  equilibrium  between  phases;  (c)  use  of  a  single- 


40 
phase  friction  factor  in  the   I'    term  and  substitution  of  mean  specific  volume 
values;  and  (d)  use  of  a  weighted  viscosity  in  the  evaluation  of  Reynolds  number 
for  determination  of  the  single-phase  friction  factor.    The  mean  viscosity  is 
usually  determined  by  the  relationship 

_i_  =  U   ~  x>    +JL.    .  (68) 

^m  «i  »g 

Although  most  investigators  would  suspect  that  this  condition  is  seldom 
fulfilled,    many  useful  results  have  been  derived.     Generally,    the  difficulties  of 
utilizing  this  method  involve  the  proper  selection  of  the  single -phase  friction 
factor,    as  well  as  the  determination  of  physical  equilibrium  between  phases. 
Investigators  presently  evaluate  the  friction  factor  from  standard  correlations 
for  liquid  flow,    using  Reynolds  number  based  upon  the  liquid  viscosity  mentioned 
above.     For  forced  convection  flow,   Owens  (40)  shows  that  the  two-phase  friction 
factor  may  be  assumed  to  be  the  same  as  that  for  the  liquid.     The  assumption 
of  physical  equilibrium  is  undesirable  but  necessary  to  effect  results. 

The  most  complete  expression  for  [Try- 1         may  be  derived  by  using  the 
momentum  equation  (Eq.  53)  and  the  following  expressions  for    a  and  da: 

a   -    (1   -X)'v/Vxv  (69) 

i  g 

v     v     dX  +  X  (1   -  X)  v     dv 
da  =     l     g  —A — ^  (70) 


[(1   -X)  v     +Xv  J 


This  derivation,    however,    is  quite  algebraically  involved.     A  more  managable 
approach  by  Owens  (40)  uses  a  mean  specific  volume  written  from  the  continuity 

equation  (assuming  the  velocities  of  the  two  phases  are  equal). 

w     v     +  w     v 

v  =      g     S i — s-  (71) 

w  x      ' 

=  Vl  +  X<Vg'vi)  (72) 
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Owens'    components  of  Eq.    41  are: 


( 


dP\    _   fTP 


G2v 


dL/, 


dDg 

(Both  Owens  (40)  and  Dukler  (14)  discuss  the  problems  of  selecting  a  proper 
two  phase  friction  factor.) 


(73) 


dP 
dL 


)      =  — 


2 
c    L 


45_dx  .  £v_  dWaPX 


dX  dL         <5vtf  dP  VOL/ 


TPJ 


(74) 


Using  the  differential  of  Eq.    72  and  substituting  into  Eq's.    (73)  and  (74), 
Owens  obtained  for  a  final  expression  for  two  phase  pressure  drop: 


■(*) 


fTP  °2  v1 

2  Dg„ 


T 
I 


1  +  Xi 


^ 


\\    G   yA  r 
1     + 


/J 


g. 


dX 


v,       X  I 


TP 


1  + 


XQ2      dvg 


(75) 


dP 


To  consider  elevation,   a  [  -j^-  \      term  must  be  added  to  the  numerator  of  the 
above  equation. 


ldL;z 


dP\ 
dLiz 


sin  0 


r 


V, 


1  +  X 


/  v 
v„ 


)] 


(76) 


Stepwise  integration  over  the  desired  interval  may  now  be  performed  to  obtain 

a  final  answer  which  will  be  dependent  on  how  X  varies  with  L,    and  the  validity 

v  dv 

of  fr^-r-,.     For  the  case  where    — ^  -  1,   and  — rl^     are  essentially  constant,   and  for 
TP  v„  dP  J 


a  linear  variation  of  X,   Eq.    7  5  becomes: 


v. 


/I  +  X2B  i    ,  AL  i  dpi  A 


/L    v  Ai 
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(77) 


ALA     dP\ 

B  H^ 


where 


vff  dv  2  2 

A  =  -^  -  1,     B  =  -nf    —  .    and  C  =  —  v„. 
^  dP      gc  g        jg 

Other   "friction  factor"  methods  (derived  by  analogy  from  single-phase 
flow  in  which  the  two-phase  pressure  drop  can  be  expressed  by  the  conventional 
Fanning  equation)  have  been  discussed  in  a  paper  by  Govier  and  Omer  (18).     In 
general,    "friction-factor"  methods  do  not  require  the  assumptions  made  in  the 
special  case  of  the  homogeneous  model,    and  can  be  applied  equally  well  to 
vertical  or  horizontal  flows;    such  methods  of  "correlation  by  analogy"  give 
good  results  even  though  the  choice  of  correlating  parameters  varies  widely. 
Excellent  papers  reviewing  this  correlation  have  been  presented  by  Owens  (40), 
Bennett  (9),   Isbin  et  al.    (27),    and  Tong  (48). 

Annular  Flow  Models 

Annular  flow  occurs  when  the  gas  travels  up  the  tube  at  a  high  velocity 
and  the  liquid  forms  an  annular  film  around  the  tube  walls.     Methods  which 
predict  pressure  drop  using  this  model  have  been  widely  discussed  due  to  the 
popularity  of  Martinelli's  et  al.    (35)  early  correlations  based  on  the  annular 
flow  model.     For  a  complete  listing  of  papers  presented,    reviews  by  Isbin  et  al. 
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(27),  Scott  (43),   Bennett  (9),   Jens  and  Leppert  (30),    Tong  (48),    and  Gresham 
et  al.    (20)  provide  excellent  summaries  and  comparisons. 

Briefly,    the  Martinelli  method  originated  from  a  series  of  studies  of 
isothermal  two-phase  pressure  drop  in  horizontal  pipes.     The  investigation  was 
confined  to  frictional  pressure  drop  alone,    as  the  use  of  horizontal  pipes 
eliminated  hydrostatic  head  terms  and  isothermal  conditions  allowed  accelera- 
tion pressure  drop  to  be  considered  negligible.     The  isothermal  pressure  drop 
correlation  was  then  extended  by  a  process  of  integration  to  two -phase  flow 
with  vaporization  by  application  to  the  forced  circulation  of  water.     A  number 
of  other  correlations  between  frictional  pressure  drop  and  properties  of  the 
system  have  been  proposed  for  isothermal  gas  liquid  flow,   but  only  the 
Martinelli  method  has  been  extended  to  two-phase  flow  with  vaporization. 
Modifications  by  Levy  (32),    Marchaterre  (34),    and  Davis  and  David  (13)  based 
on  a  momentum -exchange  model  are  summarized  in  Scott  (43). 

Lockhart-Martinelli  correlation 

This  derivation,    based  on  limiting  assumptions,   has  been  applied  to  all 
regions  of  two -phase  flow  both  by  the  originators  and  by  many  other  investi- 
gators.    Briefly,    the  original  assumptions  as  applied  to  isothermal  two-phase 
pressure  drop  in  horizontal  pipes  were:    (a)    the  static  pressure  drop  for  the 
gas  phase  is  equal  to  that  of  the  liquid  phase,    regardless  of  flow  pattern;  (b) 
the  volume  occupied  by  the  gas  plus  that  occupied  by  liquid  at  any  instant 
equals  the  pipe  volume;  and  (c)  the  frictional  gas  pressure  drop  equals  the 
frictional  liquid  pressure  drop  which  equals  the  static  pressure  drop,    i.e., 
the  momentum  and  hydrostatic  pressure  drops  are  negligible  compared  with 
the  frictional  pressure  drop;    (d)  two-phase  flow  can  be  divided  into  four  types 
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which  are  combinations  of  the  two  phases  with  either  phase  in  viscous  or  turbulent 
flow.     It  was  assumed  that  laminar  and  turbulent  flow  could  be  defined  for 
each  phase  as  the  usual  function  of  Reynolds  number  based  on  1he  particular 
phase  flowing  in  the  pipe  alone. 

Martinelli  et  al.    (35)  originally  introduced  correlations  for  turbulent- 
turbulent  and  viscous -turbulent  (liquid-gas)  flow.      Martinelli  et  al.    (37)  later 
extended  their  correlation  to  viscous -viscous  flow.     Finally,    all  four  com- 
binations,   as  well  as  for  void  fractions,    were  presented  by  Lockhart  and 
Martinelli  (33). 

Essentially,   the  correlation  presents  a  plot  of  an  empirical  function,  4>, 

against  a  parameter,    X,    with  one  curve  representing  each  of  the  four  flow 

regimes.     The  correlating  quantities  are  defined  as 

9       (AP/AL)Tp 
*g         (AP/AL)g 

or 

(AP/AL)^ 

**       (ap/au^ 

0       (AP/ALL 

2=- ' -I  (80) 

x         (AP/AL) 

where  X     represents  the  ratio  of  the  pressure  drop  functionals  of  each  phase, 
i.  e., 

«*-£)  m)  -m  m 

for  the  turbulent  situation. 

In  Eqs.    (78)  and  (79),    the  quantities  (AP/AL).  or  (AP/  AL)     are  calculated 

*■  g 

from  conventional  single-phase  correlations  on  the  basis  that  the  liquid  or  gas 


and 
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is  flowing  in  the  pipe  alone  at  the  same  individual  mass -flow  rate  as  the  two- 
phase  case.     Complete  expressions  for  X    ,    and  the  relationships  between  <J>   , 

O 

<j>^,    R       and  R     are  given  by  Lockhart  and  Martinelli  (33).     The  relationships 
are  shown  graphically  in  Fig.  4. 


^  1.00 


-*  0.10 


0.01 


Z 

o 


Fig.  4. 


1.00 
PARAMETER  X 


Lockhart  and  Martinelli' s  faired  curves  showing  relation  among 
4>o,  §„,    R*,  and  Rg  for  all  flow  mechanisms. 


Originally,    the  validity  of  the  correlation  was  questioned  for  high  and 
low  values  of  the  parameter  X,    and  for  two-phase  mixtures  near  the  critical 
point.     Later  investigations  revealed  that  the  calculation  of  pressure  drops 
by  the  Lockhart -Martinelli  method  appear  to  be  reasonably  useful  only  for 
the  turbulent -turbulent  regions.     Although  it  can  be  applied  to  all  flow  patterns, 
accuracy  of  prediction  will  be  poor  for  other  cases.     Perhaps  it  is  best  con- 
sidered as  a  partial  correlation  which  requires  modification  in  individual 
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cases  to  achieve  good  accuracy.     Certainly  there  seems  to  be  no  clear  reason 
why  there  should  be  a  simple  general  relationship  between  the  two-phase 
frictional  pressure-drop  and  fictitious  single-phase  drops.     At  the  same  value 
of  X   in  the  same  system,    it  is  possible  to  have  two  different  flow  patterns 
with  two-phase  pressure  drops  which  differ  by  over  100%.     The  Lockhart- 
Martinelli  correlation  is  a  rather  gross  smoothing  of  the  actual  relationships. 
To  illustrate  the  application  of  the  Lockhart-Martinelli  correlation,    an 

example  problem  adapted  from  Lockhart-Martinelli' s  original  paper  (36)  is 

14 
given. 

"it  is  desired  to  estimate  the  pressure  drop  per  foot  of  tube  for  the 

flow  of  oil  and  air  in  a  capillary  tube  0.00488  ft  in  diameter.     The  properties 

of  the  oil  and  air  and  the  flow  rates  are  given  below: 

w     =  2.55  X  10"5  lb/ sec 
g  ' 

p     =  0.075  lb/ft3  (at  1  atm,    70°F) 

u     =  1.22  X  10~5  lb/ sec  ft 
g 

w^  =  3.07  X  10"5  lb/  sec 

Pjl  =  54.3  lb/ ft3 

\i     -  4.25  X  10~3  lb/ sec  ft 

Let  us  assume,    for  the  moment,    that  the  flow  is  viscous -viscous,    items 

1  through  4,   below,    are  necessary  to  check  this  assumption. 

1)    Calculation  of    X      : 

w 

From  Table  2,    first  row,    first  column,    it  follows: 

2         Wl  pgMl  2.07  X  0.075  X  4.25  X  1Q2        n  „__  (R9, 

vv       w„  p.  u  2.55X54.3X1.22 

£    a     g 

.-.   \         =  0.62D 

vv 


14 

Arranged  from  the  unpublished  notes  of  Mr.    P.  S.  Martino,    Kansas  State 

University  Graduate  School. 
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and 


2)  Calculation  of  4>        ,   9.      ,    R.  and  It    : 

gvv       Jlvv        I  g 

From  Fig.  4,    with    y        =  0.620: 

<J>  =  1.98 

gw 

R„   =  0.1782  and  \\      =  0.8218 
I  g 

From  Table  2: 

^vv"    X  -0TG2U    "  3"15  (83) 

vv 

3)  Calculation  of  D„  and  D    : 

i  g 

From  Table  2,    last  column,    last  two  rows: 

D  f  2  ? 

D2"      =  ^  4>ivv    =  0.1782  (3. 15)^  =  1.77  (84) 

2 

D2       =  Rg^|vv    =  0.8218  (1.98)2  =  3.22  (85) 


Therefore 


D     =  0.00367  ft 


D     =  0.00272  ft 
g 

4)    Calculation  of  Re„  and  Re    : 

^  g 


C                     C 
f     =    S-  =  S —  (86) 

6        Rem      /4        wg       \m 

C„  C 

f*'^F 'V'i  -■■  <87> 

77  "'ue"i 
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From  Eqs.    (86)  and  (87),    with  n  =  m  =  1, 

4w 
Reg  -  TD^JjT  <86a) 

4w 
Re.  = p-p-^ —  (87a) 

But  from  Table  2, 

avv  =  R^/vv  =  01782  X  1"77  =  °-314  <88> 

Kv  =  Rg2  *|w  =  °-8218  X  3-22  =  2-64  <89) 

Thus 

Re     =  4X2-55X10"5 =  370 

g      3.14  X  2.64  X  2.72  X  10~3  X  1.22  X  10~5 

Similarly, 

Re     =  5.38. 

With  these  magnitudes  of  the  Reynolds  numbers,   assume  viscous- 
viscous  flow.     (Viscous  flow  -  R    <1000,    p.  12.) 

n 
5)    Calculation  of  the  pressure  drop: 

Once  the  viscous -viscous  character  of  the  flow  has  been  established, 

the  pressure  drop  is  readily  calculated. 

From  9 

„/4f  ~n      '      n       2-n 

(API   .    2l¥        c**iwi  (90) 

with  n  =  1  and  C  ^  =  C  '    =  16, 


8X16^    Wjg 
^L 4~  -  0.025  psi/ ft 


Jt       *  Dp    Pje  gc 
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From  Table  2 

(lr)TP  =  C  (si:),  =  (:MS>2  <002r'>  '  °'248 "*/« 

The  measured  pressure  drop  was  0.295  psi/ft,    and  the  deviation  is  of 
the  order  of  6%." 

Martinelli -Nelson  correlation 

According  to  Jens  and  Leppert  (30)  in  1955,    the  following  quotation 
provides  an  interesting  opinion: 

"It  is  the  opinion  of  the  authors  that  the  most  generally  useful  method 
available  at  the  present  time  for  predicting  pressure  drop  during  forced- 
circulation  boiling  in  tubes  is  the  one  proposed  by  Martinelli  and  Nelson  (36). 
This  statement  is  made  in  spite  of  the  fact  that  the  method  represents 
essentially  an  extrapolation;  not  of  experimental  data  for  boiling  flow,    but  of 
data  for  adiabatic  flow  of  gas -liquid  mixtures  (air  and  water,    for  example). 
Furthermore,    the  prediction  is  for  maximum  and  minimum  limits  of  pressure 
drops  which,   in  many  practical  cases,    are  separated  by  a  factor  of  two  or 
more." 

These  empirical  correlations  were  originally  based  on  data  obtained 
for  isothermal  horizontal  flow  at  pressures  close  to  atmospheric  (to  50  psi), 
normal  temperatures,    and  pipe  diameters  to  1  inch,    using  air  and  eight 
different  liquids.     Essentially,    the  paper  extends  the  Lockhart -Martinelli 
correlation  to  include  conditions  in  which  the  vapor  flow  rate  increased  and 
the  liquid  flow  rate  decreased  directly  with  flow  length,    as  in  the  forced 
circulation  boiling  of  water.     The  basic  assumptions  made  were  as  follows: 
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a)  at  point  conditions  in  the  boiling  tube,   thermodynamic  equilibrium 
exists; 

b)  the  flow  type  was  turbulent -turbulent  and  the  original  correlation 
between  <J>         and  X      could  be  used  as  a  basis  for  correlation; 

c)  on  the  basis  of  a)  and  b),    point  conditions  in  the  tube  could  be 
described  by  the  relationship 

&L  a  ♦« 

d)  a  linear  relationship  exists  between  steam  quality  and  tube  length. 
The  Lockhart-Martinelli  correlation  was  then  integrated  with  pressure 

level  as  a  parameter.     The  resulting  values  did  not  agree  with  experiments 
at  intermediate  pressures,    or  with  requirements  at  the  critical  point,    so 

corrections  were  applied  to  the  calculations.     In  effect,    an  empirical  de- 

2 
pendence  was  enforced  on  the  factor  <\>    .     Similarly,    given  values  for  R     or 

2 
R     as  a  function  of  X     at  atmospheric  pressure  and  at  critical  conditions 

(R     =  R.),    values  for  other  pressures  were  arbitrarily  interpolated.     The 

resulting  set  of  curves  given  by  Martinelli  and  Nelson  (36)  allow  calculation 

of  frictional  pressure-drops  in  the  steam-water  system  (Fig.  5). 

The  momentum  pressure  drop  can  often  be  neglected,    particulary  when 

no  mass  transfer  occurs  in  a  system.     Three  equations,    all  approximate, 

are  given  below  for  calculation  of  these  acceleration  pressure  losses,  A  Pa  , 

between  two  sections,    1  and  2,   with  the  mean  velocities  of  gas  and  liquid, 

2 


G 

a  a 


X22  X2  d-x2>2  U-Xp2 


Pg2Rg2       PglRgl       Piad-V       Pil(1"VJ 


(93) 


For  an  inlet  quality,   X     =  0,   ( Mar tinelli -Nelson  form) 


A  P. 


m 


X. 


g, 


g2*g2 


(i  -x2r 


For  an  assumption  of  non-slip  or  homogeneous  flow, 


r„G 

2     m 

g„ 
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(94) 


AP     =  — 
a       g 


m 


Xr 


g2 


X. 


gl 


+ 


*2 


X 


'j02      ^1 


(95) 


Graphs  giving  the  quantity,    r?,    as  a  function  of  quality    and  pressure  have 
been  given  by  Martinelli  and  Nelson,    based  on  their  empirical  correlations 
(see  Fig.  6). 
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Fig.  5.    Martinelli  and  Nelson's 

plot   of  c}>^tt,  <}>    u,    R^  and  Rg, 

versus  ^/x" 


Fig.  6. 


it 


100  1000  10,000 

PRESSURE-psia 


Martinelli  and  Nelson's 
multiplier  r2  versus  abso- 
lute pressure  for  various 
exit  qualities  [the  pressure 
drop  due  to  acceleration  of 
fluid  during  evaporation 
equals  r2  (C.2/g).  ] 
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For  no  mass  transfer  between  phases,    the  acceleration  loss  is 
approximately 


AP 


G 


1 


1 


pi2Ki2       PnRn 


+  G 


1 


1 


R 


R 


P    ,0    ri       r,  P        1     ft       * 

g2     g2         Hgl     gl 


(96) 


A  review  by  Sher  (44)  discusses  the  Martinelli  prediction  of  two-phase 
densities.     An  excellent  summary  of  both  Martinelli  methods  was  presented 
by  Lockhart  and  Martinelli  (33)  in  1949. 

Asa  further  aid  to  assist  in  using  the  Martinelli  method  to  predict  two- 
phase  pressure  drop,    a  sample  computer  code  (from  the  unpublished  notes 
of  G.    W.   Patraw)  is  presented  in  Appendix  B. 

Mixed  Models 

Mixed  models  are  characterized  as  having  features  based  upon  both 
the  Martinelli  and  Homogeneous  models,    even  though  the  assumptions  for  both 
are  contradictory.     Homogeneous  flow  is  implied  by  the  assumption  of  no 
radial  variation  in  velocity  and  density;  however,    the  Martinelli  correlation 
is  used  for  two-phase  pressure  drops.     In  some  cases,   the  choice  of  con- 
ditions is  based  upon  convenience,    or  is  made  in  an  attempt  to  bracket  flow 
conditions.     From  the  very  nature  of  the  premises  made,    however,    it  is 
expected  that  mixed  models  can  serve  only  a  very  limited  application. 

The  Harvey-Faust  model  (22)  applied  to  steam-water  flow,    and  one 
which  was  also  adopted  by  Rogers  (41)  to  the  flow  of  vapor -liquid  hydrogen, 
represents  an  extreme  choice  of  assumption.     When  velocity  and  density 
are  assumed  to  be  homogeneous  throughout  the  mixture  for  the  Homogeneous 
method,    by  using  the  Martinelli  correlation  for  two-phase  pressure  drop, 
the  results,    though   useful,    still  should  not  be  expected  to  replace  exact 
correlations  which  apply  to  specific  problems. 
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1  5 
Approximation  Theory       (Levy) 

Recently  an  approach  through  the  solution  of  the  equations  of  motion 
was  presented.     Using  Prandtl  mixing-length  theory  and  the  equations 
for  potential  flow  to  develop  an  equation  for  the  density  distribution,    Levy 
(32)  solved  the  equations  of  motion  by  numerical  means.     The  resulting 
relationships  between  shear  stress  and  the  operational  variables  showed 
agreement  only  with  very  restricted  data. 

Of  the  available  analytical  solutions,   the  mixing  length  model  is  most 
akin  to  the  one  proposed  by  Bankoff  (8).     Bankoff  treated  the  two-phase  flow 
as  a  continuous  medium  and  postulated  velocity  and  void  profiles  of  the  power 
type.     Levy's  (32)  solution  confirmed  Bankoff's  basic  assumption  about  the 
void  and  velocity  distribution.     Forms  of  the  mixing -length  theory  differ 
from  Bankoff's  in  that  the  mixing-length  solutions  were  found  to  vary  with 
flow  conditions  and  geometry. 

Generally  speaking,    Levy's  model  reduces  the  two-phase  system  to  a 
single -phase  system,    and  the  numerous  analytical  methods  previously 
developed  for  single -phase  flow  are  then  applied  to  solve  corresponding  two- 
phase  problems.     Because  of  the  rigorous  numerical  and  theoretical  cal- 
culations involved,   the  real  value  of  Levy's  model  lies  not  in  its  application 
to  two-phase  pressure  drop  but  in  its  applicability  to  several  other  two- 
phase  problems  which  to  date  have  escaped  solution  (i.  e  .,    entrainment  of 
liquid  by  gas  streams,    rough  channels,    entry  length  in  two-phase  flow,    heat 
transfer  in  two-phase  flow,    and  flow  with  a  gas  core  where  the  density  at  the 

center-  of  the  channel  p       falls  bo  low  the  gas  density  p    ). 

'm  &  g 


1 5 
'  Approximation  theory    involves  the  procedures  by  which  complex 

equations  may  be  simplified  to  give  an  approximate  solution 
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DISCUSSION 

A  number  of  correlation  comparisons  have  been  presented  in  the 
literature.     Various  techniques  have  been  used  to  evaluate  these  correla- 
tions.    The  most  recent  is  presented  in  a  paper  by  Dukler  et  al.    (15)  of  the 
University  of  Houston.     A  computer  recorded  experimental  data  taken  from 
published  and  unpublished  sources.     Approxiately  15,000  data  points  were 
collected.     A  coding  system  was  then  developed  to  permit  rapid  location 
covering  any  specified  range  of  conditions.     The  data  were  grouped  accord- 
ing to  test  section  characteristics,    run  conditions,    and  experimental  results. 
The  resulting  deck,    designated  the  "data  bank,  "  is  used  to  compare  the  data 
resulting  from  relatively  consistent  experimental  observation  with  the  results 
calculated  from  the  various  methods  available.     By  statistical  calculations, 
each  method  may  then  be  compared  for  accuracy  by  utilizing  three  statistical 
parameters:    d,   the  average  percent  deviation  of  the  calculated  values  from 
the  measured  values;  a,   the  variance  of  individual  percent  deviations;  and  % 
which  is  the  width  of  the  band  around  d  that  contains  68%  of  the  calculated 
d' s.     Five  correlations  were  compared  by  Dukler,    who  used  line  size, 
viscosity,   pressure     range,    and  flow  pattern  as  variables  (Tables  3,  4,  5). 
The  overall  comparison  indicated  the  Lockhart-Martinelli  correlation  was  better 
than  the  other  four  except  for  two  special  cases.    Regression  analysis,    in- 
cluding a  comparison  of  scatters  between  Martinelli's  correlation  and  the 
experimental  data,    verified  that  the  method  needs  still  more  refinement. 

Most  investigators  have  chosen  to  compare  their  methods  to  the  original 
Lockhart-Martinelli  correlation  or  the  Homogeneous  method.     Reviews  by 
Isbin  et  al.    (28)  and  Gresham  et  al.  (20)  compare  each  correlations' s 
theoretical  and  physical  limits.      Anderson  and  Russell  (3)    using 
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the  same  method  as  Dukler  (Tables  6,  7),    furnish  additional  correlation 
comparisons.    Hughmark's  comparisons  (2fi)  are  presented  in  Table  8. 

Brief  mention  should  be  made  concerning  the  present  statvis  of  inves- 
tigation.    Journals  published  monthly  by  ASME  and  AIChE  provide  good 
sources  of  information  and  periodically  list  special  conferences  held  to  dis- 
cuss the  status  of  current  or  recent  projects.     A  review  of  the  proceedings 
are  frequently  published  in  professional  journals.     Since  improvements  in 
two-phase  systems  depend  considerably  upon  break-throughs  in  other  areas, 
studies  relating  to  mass  transfer,    heat  transfer,    chemical  reaction,    flow 
stability,    critical  flow,    and  flow  over  restrictions  are  especially  important 
to  investigators  working  on  new  pressure-drop  correlations.     Complete 
discussions  of  these  areas  may  be  found  in  Tong  (48),    Scott  (43),    Anderson 
and  Russell  (3),    and  Isbin  et  al.    (28). 
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Table  3.     Dukler's  comparisons  of  pressure  drop  correlations  for 
various  pressure  ranges  (one -component  system). 

Pressure  Clicnoweth  l.ockliart 

Hanne  D  Maker        Mankol'f  Martin Marunolli  Ya^fi 

(psia)        (in.)      d        or      *        d  a  *  d  a  v  u  a  &  d  a 


25- 
100 

400- 
800 

1000- 
1400 

All  Data 
Points 


1        234  230  -  4,852  3,904  -  53.5  05.9  07.5  39.5  66.6  55.0  -38.8  65.2  45.0*  93 

1/2     632  368  -  377  206  -  55.9  88.2  65.0  63.6  57.0  52.5*  -87.7  8.8       7.5  103 

1/2     358  177  -  31.9  73.0  52.5  63. G  145  112  -17.5  21.1  17.5*  -05.6  16.2  10.0  131 

414  310  6.6  128  23.5  GO. 7         *  -72.6  42.6  327 


*  Correlation  seemed  to  predict  best  for  each  row. 


Table  4.     Dukler's  comparisons  of  pressure  drop  correlations'   showing 

effects  of  line  size  and  liquid  viscosity  (two-component  systems). 


CI 

lenow  e 

Hi 

1 

liart 

D 
(in.) 

(cp) 

Maker 

Mankol'f 
d          a       9 

Martin 

A 

la  il  iiu 

•1H 

Yaci 

n 

CTD 

d 

a 

* 

d 

<r 

>!' 

d 

a 

* 

d 

a 

* 

D 

1 

64.2 

40.0 

45.0 

2000 

900 

- 

-8.5 

17.0 

1  5.0 

-  6.  6 

10.1 

1  0.  0* 

40.9 

29.1 

30.0 

224 

6.0 

5.0 

1 

3 

77.4 

335 

87.5 

1172 

2220 

- 

11.2 

55.6 

30.0 

3.8 

29  '. 

2  0.0  * 

103 

123 

- 

230 

22.0 

17.0 

20 

30.7 

09.5 

40.0 

737 

1384 

- 

42.5 

94.2 

65.0 

-5.5 

24.7 

2  0.0* 

401 

192 

- 

156 

18.0 

14.0 

1 

-13.6 

60.3 

65.0 

1170 

2910 

- 

-2.7 

24.8 

2  0.0* 

9.2 

37.7 

2  5.0 

1G2 

228 

- 

320 

16.0 

7.0 

2 

3 

19.3 

79.0 

82.5 

4810 

4654 

- 

8.4 

45.3 

45.0 

-4.7 

22.!' 

25.0* 

62.3 

74.5 

80.0 

398 

16.0 

9.0 

20 

73.0 

159 

90.0 

2004 

4893 

- 

95.4 

268 

- 

13.2 

52.9 

30.0* 

271 

32  5 

- 

401 

24.0 

8.0 

1 

11.5 

79.2 

82.5 

21  7G 

3072 

- 

15.0 

4  0.2 

30.0^ 

31.0 

50.2 

47.5 

27.6 

104 

97.5 

109 

- 

- 

3-1/2 

3 

7.1 

60.1 

72.5 

4720 

5000 

- 

27.0 

G2.0 

45.0 

1  6.  3 

3  9.3 

22.5* 

84.5 

86.1 

- 

67 

18.0 

12.0 

20 

31.8 

60.6 

47.5 

2432 

3561 

- 

51.0 

91.0 

57.5 

-0.4 

2  C).  2 

22.5* 

147 

83.4 

- 

111 

27.0 

13.0 

1 

-70.5 

11.6 

10.0 

254 

213 

- 

51.2 

23.7 

30.0 

30.3 

12.2 

12.5* 

93.3 

22.3 

22.5 

24 

- 

- 

5-1/2 

3 

-0.5 

44.6 

45.0 

209G 

3704 

- 

20.0 

57.5 

4  5.0 

11.6 

-1 1 . 5 

37.5* 

10G 

80.8 

- 

131 

19.0 

13.0 

20 

7.G 

47.8 

50.0 

2692 

5263 

- 

37.1 

79.4 

47.5 

-1.0 

2  4.0 

2  5. 0* 

120 

69.0 

- 

122 

21.0 

7.0 

Ail 

Data 

28.2 

159 

G5.0 

27.  6 

0  9.5 

42.5 

4.2 

5G.0 

2  5.0* 

155 

104 

2293 

Points 

T  d,  a, 

and 

tt  expressed  as  per 

centat 

;es. 

*  Correlation  seemed  to  predict  best  for  eacli  row. 


58 


Table  5.     Dukler's  comparison  of  pressure  drop  corrections  showing 
effects  of  flow  pat  lorn. 


Observed  Chenoweth  Lockhart 

flow  Baker  Martin  Martinelli 


pattern  da  *  n         d  a  #  d  a  * 


nt 


Bubble               15.6  168          30.0  960       -           -            -  -  -           -  0 

Plug  116  92.9  100  69      -6.5      19.8  10.0*      9.4  36.3   20.0  270 

Stratified           -           -            -  -      122        103           -  23.3  33.0   22.5*  34 

Wave  -91.0        2.5        2.5  98    138        301  135  38.4  85.7    42.5*  287 

Slug                     61.0  218  135  251      12.8      92.1  27.5  2.9  31.2    17.5*  974 

Annular             68.7  81.2  100  430      22.2      81.1  65.0    -12.8  35.6    30.0*  265 

Dispersed        16.9  35.0      30.0  133      84.4      81.1  85.0  18.0  34.1    25.0*  111 

t  Except  for  Baker's  Correlation,   which  is  based  on  the  number  of  points 
calculated  (14)  to  be  in  each  flow  regime,    all  correlations  are  tested  on  the 
same  points. 

*  Correlation  seemed  to  predict  best  for  each  row. 

Table  6.     Anderson  and  Russell's  values  of  d,   a,    and  \&  f or  pressure  drop 
correlations  for  different  line  sizes  and  liquid  viscosities 
(using  Dukler's  method). 

Line      Liquid  Dukler  et  al.  Lockhart-  Chenoweth- 

Size    Viscosity  (general) Martinelli Martin  Baker 

ln-             cP. g           g          <it         d          q          -tj  d         q          g  d            a  ¥ 

1            -25.2    10.2  13.5  -6.6  10.1  10.0  -8.5  17.8    15.0  64.2    40.0  45.0 

1  3                8.6   24.8  12.0  3.8  29.1  20.0  11.2  55.6    30.0  77.4   335  87.5 
20                6.7    24.4  18.6  -5.5  24.7  20.0  42.5  94.2    65.0  30.7    89.5  40.0 

1  2.4    18.4    15.5      9.2    37.7    25.0    -2.7    248     20.0    -13.6    60.3    65.0 

2  3  1.6    19.7    16.0    -4.7    22.9   25.0      8.4   45.3    45.0      19.3    79.0    82.5 
20              10.3    27.2    20.0    13.2    52.9    30.0    95.4   268        -         73.0    159     90.0 

1  -0.3  26.8  26.2    31.0    50.2  47.5  15.0  40.2  30.0  11.5  79.2  82.5 

3-1/2           3  9.3  24.9  25.0    16.3    39.3  22.5  27.8  62.0  45.0  7.1  60.1  72.5 

20  10.6  24.5  18.6    -0.4   26.2  22.5  51.0  91.8  57.5  31.8  60.6  47.5 

1  50.6  18.8  19.3    38.3    12.2  12.5  51.2  23.7  30.0  -70.5  11.6  10.0 

5-1/2           3  11.2  19.2  16.0    11.6    41.5  37.5  20.0  57.5  45.0  -0.5  44.6  45.0 

5*               20  7.3  22.7  14.0    -1.0   24.8  25.0  37.1  79.4  47.5  9.6  47.8  50.0 


*Size  used  by  Dukler  et  al. 
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Table  7.     Anderson  and  Russell's  values  of  d,   a,    and  ^  for  pressure  drop 
correlations  for  different  flow  patterns  (using  Dukler's  method). 


Dukler  et 
(genera 

al. 
1) 

Lockhart- 
Martinelli 

Chenoweth- 
Martin 

Baker 

Flow  type 

d 

O" 

^ 

d 

o          * 

d 

a 

^ 

d 

o 

<k 

Plug 

9.5 

18.0 

15.1 

9.4 

36.3   20.0 

-6.5 

19.8 

10.0 

116 

92.9 

100 

Stratified 

13.4 

30.3 

17.5 

23.3 

33.0   22.5 

122 

103 

— 

— 

— 

— 

Wave 

11.5 

22.6 

20.5 

33.4 

85.7    42.5 

138 

301 

135 

-91.0 

2.5 

2.5 

Slug 

9.5 

21.7 

15.6 

2.9 

31.2    17.5 

12.8 

92.2 

27.5 

61.0 

218 

135 

Annular 

-11.2 

26.0 

19.0 

-12.8 

35.6    30.0 

222 

81.1 

65.0 

68.7 

81.2 

100 

Dispersed 

14.8 

16.9 

17.1 

18.0 

34.1    35.0 

84.4 

81.1 

85.0 

16.9 

35.0 

30.0 
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Table  8.     Hughm  ark's  comparison  of  experimental  horizontal  flow  data 
with  correlations. 


Average  absolute  deviation, 

-\rTP,  % 

No.    of  tests 

Hugh-mark'  s 

Lockhart- 

Chenoweth- 

16 

for  Hughmark 

series  no. 

Martinelli 

Martin 

Pocttm; 

an     Hughmark 

method 

1 

19.2 

14.9 

.;;; 

15.5 

59 

2 

38.6 

1<).  9 

32.0 

14.3 

77 

3 

21.9 

23.2 

72 

21.6 

61 

4 

75.3 

18.5 

35.2 

24.3 

46 

5 

18.6 

20.1 

30 

26.4 

22 

6 

21.7 

15.6 

45 

28.8 

33 

7 

19.4 

— 

19.3 

25.3 

59 

8 

18.7 

— 

22.3 

12.3 

33 

9 

31 

— 

25.5 

12.2 

20 

10 

17.2 

22.3 

43 

26.7 

50 

11 

18.2 

25.2 

20.7 

20.8 

18 

12 

29.6 

— 

14.3 

27.9 

22 

13 

25.7 

20.8 

17.5 

21.6 

28 

14 

28 

26.7 

55 

27.5 

25 

15 

2  4.1 

49 

83 

18.9 

228 

16 

13.0 

6.3 

30.5 

11.7 

15 

17 

19.7 

8.0 

25.5 

6.0 

18 

18 

2  5.2 

4.8 

23.5 

9.2 

18 

19 

42.3 

19.6 

49.5 

24.5 

29 

20 

22.7 

22.6 

27.4 

10.1 

40 

21 

32.6 

— 

46.5 

16.9 

31 

22 

90.5 

43 

55 

42.5 

15 

23 

63 

45 

54 

38 

8 

24 

8.2 

12.7 

35.2 

6.9 

57 

25 

31.3 

37.8 

35.4 

24 

30 

26 

28.4 

34.8 

14.4 

24.6 

25 

27 

18.5 

18.1 

85.5 

23 

23 

28 

22.9 

17.8 

69.5 

18.2 

26 

29 

2  4.7 

22.4 

47 

15.3 

27 

30 

60.2 

2  6.6 

32.1 

19.5 

48 

31 

63.1 

20.1 

24.9 

13.3 

31 

32 

92 

24.7 

26.6 

11.4 

33 

33 

45 

15.3 

24 

11.6 

52 

34 

61.6 

20.8 

35.2 

10.6 

30 

35 

69 

15.2 

22.9 

8.3 

20 

36 

67 
.       36.3 

2  3.4 
22.4 

45.5 
38.4 

11.2 
18.9 

22 

Av 

Total    1379 

16Poettman,    F.   II.,   et  al.,   J.    Petrol.    Technol.    8:17  (1956) 
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SUMMARY 

This  paper  presents  a  literature  review  of  selected  two-phase 
pressure  drop  correlations,    with  enough  background  material  to  indicate  the 
major  source  of  difficulties  being  encountered  currently  in  such  an  analysis. 
From  this  review,    the  conclusions  are: 

1.  Two-phase  flow,  as  classified  according  to  system,  appearance  and 
type,  involves  so  many  variables  that  the  possibility  for  obtaining  a  universal 
solution  to  the   problem  appears  realistically  unfeasible. 

2.  Two-phase  pressure  drop  consists  of  the  summation  of  pressure 
losses  by  friction,    momentum  exchange  between  phases,   and  elevation 
phenomena.     The  components  due  to  acceleration  and  elevation  are  much 
easier  to  calculate  than  the  components  due  to  friction.     Energy,    momentum, 
and  continuity  equations  are  fundamental  to  the  analysis  of  all  three.     This 
paper  reviews  these  equations  and  includes  the  physical  significance  of  vari- 
ous simplifying  assumptions  which  permit  mathematical  solution. 

3.  Correlations  and  example  methods  have  been  presented  which  may 
be  categorized  into  the  following  four  approaches:  empirical,    similarity, 
mathematical  analysis  of  the  simplified  physical  models,   and  approximation 
theory.     The  mathematical  analysis  of  simplified  physical  models  in  the  forms 
of  the  Homogeneous  and  Martinelli  methods  have  received  greatest  attention 
and  application  of  the   newer  approaches,   Dukler's  method  appears  most 
promising.     Depending  on  the  specifics  of  the  problem,    deviation  from  experi- 
mental data  (20%  to  80%)  must  be  expected.     The  major  inability  of  the 
simplified  physical  model  approach  to  correlate  with  experimental  data  may 
be  traced  to  its  neglect  of  important  variables  and,    in  particular,    to  its 
failure  to  deal  with  the  flow  structure.     The  hydrodynamics  of  two-phase  flow 
will  require  a  great  deal  more  study  before  more  exact  results  may  be  expected. 
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RECOMMENDATIONS 

It  is  suggested  that  this  paper  be  utilized  for  preliminary  analysis  of 
two-phase  pressure  drop  by  investigators  who  are  unfamiliar  with  the  variety 
and  limitations  of  the  correlations  presented  in    the  literature. 

Pertinent  recommendations  of  the  author  are: 

1.  Greater  emphasis  should  bo  placed  on  the  evaluation  of  correlation 
methods  currently  found  in  the  literature.      It  is  now  the  responsibility  of  the 
investigator  to  reevaluate  these  correlations,    analyze  their  errors  and/ or 
limitations,    and  construct  better  correlations  therefrom. 

2.  Though  short-term  research  analyzing  methods  and  modification  of 
existing  methods  certainly  will  be  very  useful,    emphasis  should  be  placed 
upon  longer  range  programs  which  will,   in  time,    significantly  diminish  the 
inherent  complexities  and  experimental  difficulties  presently  encountered. 
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Baker's  Flow-Prediction 
Example  Problem 

This  example  will  demonstrate  the  use  of  the  flow -prediction  chart 
by  Raker. 

Given:    A  1 -in.    horizontal  pipe,    a  water-flow  rate  of  2180  lb/hr,   an 
air-flow  rate  of  10.9  lb/hr,    and  flow  conditions  such  that  the  average  physical 
properties  of  the  air  and  water  are  estimated  to  be  as  follows:    /j    =  1  cP, 
p    =  62.3  lb/ ft3,  y    =  73  dynes/cm,    pa  =  0.075  lb/ ft3. 


Find  the  expected  flow  pattern  by  using  Plate  I  and  the  method  of 


Baker: 


Calculation  of  X' 


1/2 


X1  =     (p  70.075)  W62.3) 


X'  =   [(0.075/0.075)  (62.3/62.3)  ^2 


(97) 


X'  =  1.0 


Calculation  of  i//' 


'/•' 


1/3 


73/7LJ  [^(62.3/pL) 
^'  =   [73/73)  [1(62. 3/62. 3)]1/3 


(98) 


<//*  =  1.0 


Calculation  of  G 


Gi=wiAp 


G^=  2180/  (0.25  ttX  0.0833)' 
G^  =  400,000  lb/hr-ft2 


(10) 
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Calculation  of  G   : 
g 


G     =  w  A  (9) 

g        g     P 

G     =  10.9/(0.2  5  ttX  0.0833)2 
g 

G     =  2000    lb/hr-ft2 
g 


Calculation  of  GV  <//'  /G 
J- 


G  X'tf/'/G     =  (400,000  X    1    X   l)/2,000 


G   X'^'/G     =  200 
i  '     g 


Calculation  of  G   /X'  : 

G  /X1    =  2000/  1 


G  /X'    =  2000 


From  Plate  I,   the  pattern  is  slug  flow. 
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TONG'S  SUMMARY  OF  RECENT  EXPERIMENTAL   STUDIES 
FOR  PRESSURE  DROP  IN  STEAM-WATER  MIXTURES 


Data  Range 


* 


Pressure    Mass  velocity     Quality 


Source'1'  Description  of  test  (psia)        (106  lb/hr  ft2)    (%bywt) 

Weiss  (W 13)      Electrically  heated  ver-  20-1400  0.18-1.8  0-100 

(1952)  tical  stainless  steel  tube 

0.174  in.   I.D.  X  24  in. 

long 

Isbin  et  al.  Unheated  horizontal  tubes  25-1415  0.30-1.40  3-98 

(16)  (1958)  0.484-1.062  in.   I.D.,   3ft 

8|  in.  to  8  ft  long 

Lester  (L12)      Unheated  horizontal  steel  30-100  0.04-0.54         10-85 

(1958)  pipe  lines  4.06  and  6.06  in. 

I.D.  X  40  ft 

Schrock  and       Electrically  heated  ver-  42-505  0.18-3.28  5-57 

Grossman  tical  stainless  steel 

(S7)  (1959)  tubes  0.118  in.   I.D.  X  15, 

30,   and  40  in.  long,  0.237 

in.   I.D.  X  15  and  30  in. 

long,  0.432  in.   I.D.  X  30 

in.   long 

Becker  et  al.     Electrically  heated  ver-  85-600  0.29-3.80  0.80 

(B17)  (1962)       tical  stainless  steel  tubes 

0.305  in.   I.D.   and  0.391 

in.   I.D.  X    10.25  ft 

Perroud  et  al.  Electrically  heated  and  un-    140-940  0.24-6.05  0-100 

(P3)  (1960)  heated  stainless  steel  ver- 

tical tubes  0.197  in.  I.D. 
and  0.394  in.  I.D.  ,  15.75 
and  31.5  in.   long,   respectively 

Hoglund  et  al.   Electrically  heated  ver-  150-600  0.6  0.9-6.5 

(H19)  (1961)  tical  stainless  steel  rec- 
tangular channel  2X  0.25 
X  60  in. 

Haywood  et         IX  l|  in.  I.  D.  tubes  heated    250-3000  0.54-1.26  0-57 

al.   (H8)(1961)  and  unheated  oriented  both 

horizontally  and  vertically: 

heated  lengths  16  and  24  ft 
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TONG'S  SUMMARY  OK  RECENT  EXPERIMENTAL  STUDIES 
KOR  PRESSURE  DROP  IN  STEAM -WATER  MIXTURES 


Data  Range 


Pressure    Mass  velocity      Quality- 
Source*  Description  of  test  (psia)  (106lb/hrft   )    (7obywt) 

Davidson  et        Straight  or  flat  spirally  290-3300  0.07-3.89  0-100 

al.    (Dl)(1943)  coiled  tubes;  inside  diam- 
eters from  0.5-1.75  in.  I.D., 
10  to  15  ft  long;  all  sections 
were  heated 

Cicchitti  et        Electrically  heated  and  un-    510-770  1.60-2.00        24-70 

al.  (C12)  heated  vertical  stainless 

(1960)  steel  tube  0.202  in.   I.  D.  X 
15.8  in.   long 

Armand  and       Horizontal  stainless  steel       5  10-1280  0-90 

Tretchev  (A9)  tubes,   heated  and  unheated, 
(1947)  2.2  in.   I.  D.  X  13.1  ft 

Bertoletti  et      Electrically  heated  and  un-    600-1200  0.74-2.95         20-85 

al.    (B36)  heated  stainless  steel  ver- 

(1961)  tical  tubes  and  annuli 

Tubes 

0.204  in.  I.D.X  3.9,  8.3, 
16.2,  and  31.4  in.  long 
0.126,  0.158,  0.248,  and 
0.323  in.  I.D.X  3  1.5  in.  long. 

Annuli 

0.394in.  O.D.X  0.197  in. 

I.D.X  23.2  in. 

0.394  in.  O.D.X  0.152  in. 

I.D.X  23.2  in. 

Sher  et  al.  Vertical  heated  rod  cluster  615  1.0-2.5  0-5 

(S13M1962)         assembly  16  rods  0.416O.D. 
X36in.  long  4X4  array  on 
0.570  in.  centers  (lower  half) 
0.367in.  O.D.X  36  in.  long  on 
0.570  in.  centers  (top  half)  in- 
cluding spacers  in  box  2.3  99 
in.  square 
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TONG'S  SUMMARY  OF  RECENT  EXPERIMENTAL  STUDIES 
FOR  PRESSURE  DROP  IN  STEAM- WATER  MIXTURES 


Data  Range 


Source 


* 


Description  of  test 


Pressure    Mass  velocity     Quality 
(psia)  (106ib/hrft2)    (%bywt) 


Janssen  and 
Kervinen  (J7) 
(1961) 


Columbia 

University 

Engineering 

Research 

Laboratories 

(C21)(1961- 

62) 

Quinn  (Q4) 
(1961) 


Sher  (S12) 
(1957) 


Vertical  unheated  stainless  1000  0.1-1.0 

steel  rod  cluster  assembly 
16  rods  0.423  in.  O.D.X40 
in.  long  4X4  array  on  0.76 
in.  centers  in  box  2.8X  2.8  in. 
(and  including  various  ori- 
fice fittings,  wire  spacer, 
etc .  ) 

Various  vertical  7-rod  and  1000  0.5-1.6 

19-rod  cluster  assemblies 
both  with  and  without  wire 
wraps 


Vertical  unheated  rod  cluster     1000  0.7-1.2 

assembly  25  rods  0.363  in.  O.D. 
X  29.7  in.  long  5X  5  array  on 
0.615  in.   centers  in  box  3.04X3.04 
in.    (both  with  and  without  wire 
spiral  spacers) 

Electrically  heated  stainless     1100-2000      0.7-5.0 
steel  vertical  rectangular  chan- 
nels 1.0x0.097x27  in.   and  1.0X 
0.050X27  in. 


2-26 


0-40 


0-15 


0-40 


# 


References  listed  in  the  column  are  from  Tong  (48). 
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APPENDIX  (' 
TWO-  PHASE  PRESSi   RE  DROP  COMPUTER  CODE 
(Homogeneous  and  Lockhart-Martinelli -Nelson  Methods) 

High-speed  digital  computers,^  reducing  computation  time,  now  allow 
the  engineer  to  solve  complex  design  problems.     Among  the  numerous  types 
of  problems  which  can  now  be  solved  efficiently  are  trial  and  error,   or  iter- 
ative optimization  designs.     Problems  of  this  type  normally  are  approached 
through  the  use  of  a  general  problem  code  that  relies  upon  and  calls  for 
specialized  subroutines.     These  subroutines  provide  answers  to  problems 
whose  methods  of  approach  are  likely  to  be  continually  updated  or  varied. 
Calculation  methods  of  this  type  are  especially  flexible  in  that, while  better 
subroutines  are  being  developed,  the  general  program  remains  operational. 

According  to  its  importance,   numerous  two-phase  pressure  drop 
subroutines  have  been  written  to  handle  a  variety  of  flow  conditions.     As  an 
example  of  how  the  various  computation  methods  have  been  programmed,   the 
Homogeneous  and  Lockhart-Martinelli -Nelson  subroutines,  taken  from  the  per- 
sonal notes  of  G.  W.   Patraw,   are  presented.     Tong  (48)  provided  the  prin- 
cipal reference  for  these  methods.     Both  of  the  codes  have  been  written  for 
a  6600/3600  CDC  computer  in  FORTRAN.    The  author  considers  only  the  method 
of  approach  important.     No  attempt  has  been  made  to  develop  or  refine 
these  subroutines  more  thoroughly. 
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Homogeneous  Subroutine 


»  ■     COLI 

Jt US   _J 


« 

SI  .  INE    ■      I 

L.  /DROP/    D,Y,  WiP     ,.,-!'   ,NI   Sh 

CON  /PROP/    CPrLf-    SI    SVI  SL  ,V1SG,VP,VL 

COf  :  ON/J OM/DELP 
'•■■•■  i    .     (  M  =  3.  141  5Q  )  .  f  7J  1  t  PC VI  =21  2  0.  ) 


A    =    PI*D»D/A. 

Al     =    .092*(A.*W/PI)**l;     *VISL*».2/GE/D**4.8 

A  2    =     (W/A)**2/GE 

- 

DXIDL    =    ON/Y 

-:  lsh 


DL  =Y  /AMESH 

XI    =    XIO-. 5*0X1 DL»UL 

P    =    PO*PCVT 

DO    9  9    1=1, MESH 

XI  =  XI+DXIDL*DL 

XII  •-=  1  .-XI 


VG  -  V?/P 

DVi"5  =  -VG/P 

T  =  rP+CONP/LOGF ( C P*PC .  /?) 

KEG  =  4.»XI*W/(PI*I  ■ V ! 

L  =  4.*XI1*W/(P1*I  »  VI  SL  ) 


C1C    =    A2*VG*DXlOL 
TOP    =    DPFDL-C10 
C8    =    XI*A2*DVDP 
BCT    =     1.+C8 
DP    =TOP/BOT*DL 
P    =    P+DP 


U    =    W/A*{ VL+XI*VG) 

DELPF  =  DPf DL»DL 

DELPM  =   DP-DELPF 

WRITE  OUTPUT  TAPE  ■■ ,  LO  L  ,  XI  ,  VG  ,  T  ,REG,REl  tDPFDL  ,C  10  ,  C8  ,  DPf  UiOELPF, 

X 
1001  i  '     r  (  7 E  1  2  .  1  / 5 r.  1  ?  .  '; ) 


IF     (P)     999,999,99 
99    CONTINUE 

DELP    -     l.-P/( P0*PCV1 ) 

G(  ir  555 
999  i)f  LP  =  0. 
555    R  ! 

E  . 


Lockhart-Martinelli -Nelson  Subroutine  „n 

*  CARDS  COLU 

« Li-SJ  .. 

*  r    .   -  .  ... 
SUBROUT  INE  LMN 

CI     •      */    DtYtWf  PO,  XI  ,    ,   SI 
C       /PROP/  CP,CONP, TP, VI  ....    ■  •  •  .'I- 
COMMOU/LOCK/DCLP 
DATA  (PIs3,14159).((;g«3  .    '.  .) 

a  =  pi»n*o/A. 

A  1  =  .  09 2»  (  4  •  »W/P  I  )  *»  i  .  -/:-  =  ././  -  -  • . 

a?  =  :  .:/ \  >**2/Gfc 

A3  =  SQRTF( VL )• ( VI SL/VISG )*«. 1 

ON  =  '■ 

OX  I  PL  =  Q!\/Y 

A.v  Sh 

DL=Y/AMESH 

XI  =    XIC- .5»DXIDL»I  L 
P    =    P0*PCV1 

do  99   i  =  i , ;. r  sh 

XI     =    Xl+nxiDL«PL . 

XII  =     I .-XI 

VG  =   vp/p 

DVDP    =    -VG/P 

XT    =    (  XI1/XI  )  **.9/SQK  rr-(vG)  »/ 

RXT    =    SQRTF(XT) 
Phi     =     LXPF(  1.  16«R/T  )  /X  T«      (_.      1  H»RXT)  +_XJ 

AO    =    XI 1«VL+XI*VG 

7     =    TP+CONP/LOGF(CP«PCVT/P ) 

DPFDL     =    A1*VG»XI»* 1 . 3 

DPFCL     =-PHl»*2*DPFDL 

CL            =    XI1**2»VL/ALF1 
C2    =    XI««2»VG/ALF 

C3    =    -l./XIl+DALDX*[»XrDX/ALFl 

CA    =    l./XI-DALDX»DXTDX/ALF 

THP    =    DPFDL-A2*DXIDL»(C1*C3+C2*C4) 

C5    =    C1/ALF-1-C2/ALF 

C6    =    C5*DALDX»DXTDV 
0  7    =    XI»XI/ALF 

C8    =    A2*L)VDP* (C6+C  7) 

BGT    =    1.+C8 

UP    =TOP/BOT* 

P    =    P-t-DP 

UG    =    XI *VG»W/ (ALF*A) 
UL    =     X  1  1*VL»W/  (  AL  1    1  «  A  J 

Di-LPF     =    DPFDL»DL 

DELP.'*    =       DP-DELPF 

REG    =    *r.*Xl*W/(Pl*D»VISC) 

R  E  L     =    4.*XIl«Vs/(PI*D*VlSL) 

F     =     .(K6/REG»*.2 
WRIT!      OUTPUT     T.APE      •■  ,  '  :..      .■      ,        [_,  ,     ^     .  (.PXTDV.T. 


X     CPFPL,C1,C2,C3,  C4,  C5,C6,C7,<    UPPtUGtLL.P     L        -  •         .,RtLtF 

100  1    FORMA! (7F12.     ,  .     tl     12. 3/<      I     .3) 

If"     (  I' )     999,999f99 

9?    O  ! .  I  I  Lilc 

D  L I P     =     1 .-P / ( P 

or    n    '>1^ 

999  DELP  =  C. 
555  RE1 
END 
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APPENDIX  D 

NOMENCLATURE 

2 
A  Area,    L 

C  Constant  defined  in  Blasius  Equation  (33) 

C  Wetted  perimeter,    L 

D  Pipe  inner  diameter,    L 

d  Differential 

d  Average  percent  deviation  of  the  calculated  values  from  measured 

values. 

E  Mechanical  energy  dissipated,    LP1/ M 

f  Eriction  factor 

f*         Chenoweth  and  Martin  fictitious  friction  factor  based  on  Fanning 
equation 

F  Function 

Fr        Froude  number 

2 
g  Acceleration  due  to  gravity,    L/ T 

g  32.17    LM/FT2 

&c 

G  Mass  velocity,    M/L   T 

h  Specific  enthalpy,    LF/ M 

H  Enthalpy,    LF 

J  778  ft-lbf/B 

K  Hughmark  flow  parameter.     See  Fig.  3. 

K'  Chenoweth  and  Martin  allowance  for  valves  and  fittings.    See  Eq.  55. 

K"  Bankoff's  flow  parameter.     See  Eq.  24. 

KE  Kinetic  energy,    LF 

L  Length,    L 

2 
P  Pressure,    F/ L 

3 
Q  Volumetric  flow  rate,    L   /T 
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Q'         Total  heat  flow,    LF/ T 

q  Specific  heat  flow,    LF/L 

3 

r  Martinelli -Nelson  multiplier,    L   /  M,    (36) 

R  Local  volume  fraction  (holdup) 

R  Average  volume  fraction 

Re  Reynolds  number 

5  Dimensionless  distance  from  the  wall 
u  Velocity,    L/ T 

u'  Relative  velocity,    L/ T 

U  Thermodynamic  internal  energy,    LF/ M 

3 
v  Specific  volume,    L   /  M 

_  3 

v  Mean  specific  volume,    L  /  M 

w  Mass  flow  rate,    M/ T 

We  Weber  number 

X  Quality  (vapor  mass  fraction) 

z  Elevation  distance,    L 

Z  Hughmark  dimensionless  correlating  factor,    see  Eq.  29 

Subscripts  and  Superscripts 

a  Acceleration  term 

f  Friction 

g  Gas 

6  Liquid 

m  Mixture 

m         Exponent  of  Reynolds  modulus  in  the  B  la  si  us  expression  for  the  friction 
factor  for  the  gas  phase. 

n  Exponent  of  Reynolds  modulus  in  the  Blasius  expression  for  the  friction 

factor  in  the  liquid  phasi  . 

NS        No  slip 
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o  Single  phase  flow 

p  Pipe 

s  Superficial 

t  Turbulent 

tt  Turbulent  liquid,   turbulent  gas  flow 

TP  Two  phase 

v  Viscous 

w  Viscous  liquid,   viscous  gas  flow 

W  Wall 

z  Elevation 

1,2  Position  1,    position  2 

*  Fictitious 


Greek  Letters 

a  Void  fraction 

ay  Flow  modulus  for  gas  (L-M  model) 

|3  Dimensionless  group  defined  by  Eq.    59 

/3'  Flow  type  modulus  for  liquid  (L-M  model) 

2 
7  Shearing  stress  -  F/L 

A  Small  change 

X  Ratio  of  the  volumetric  flow  rate  of  liquid  to  the  total  volumetric  flow 
rate  at  average  pressure. 

XI  Baker  flow  parameter.     See  Eq.    97. 
jj.           Viscosity,   F/TL 

it  3.1416 
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Q  Angle  of  elevation 

p  Density,    M/L 

a  Variance  of  individual  percent  deviation 

</>  Function  of  \  utilized  in  calculating  two-phase  pressure  drop 

X  Dimensionless  parameter  (Lockhart-Martinelli-Nelson  models) 

V  Width  of  band  around  d  that  contains  68%  of  the  calculated  d's 

ii  Dimensionless  group 

(//*  Chenoweth  and  Martin  fictitious  dimensionless  group 

ip1  Baker's  dimensionless  group.     See  Eq.  98. 
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ABSTRACT 

Methods  for  predicting  pressure  drops  in  two-phase  flow  systems  have 
been  investigated.     A  literature  survey  covers  the  historical  development  and 
current  status  of  the  subject.     Flow  classification,    technical  terms,    important 
equations,    and  inter  relationships    among     variables  pertinent  to  the  study  of 
two-phase  flow  have  been  reviewed.     A  general  description  of  two-phase 
pressure  drop  has  been  presented  including    the  energy  and  momentum  equa- 
tions.    Prediction  methods  have  been  categorized  according  to  method  of 
approach  (empirical,    dimensional,    similarity,    mathematical  analysis  applied 
to  simplified  physical  models,    and  approximation  theory).     Numerical  examples 
have  been  included  to  demonstrate  the  Dukler  and  Lockhart-Martinelli 
approaches.     Comparison  methods  by  Dukler,   Anderson  and  Russell,    and 
Hughmark   with  tables  have  been  included.     Baker's  chart  and  method  of  pre- 
dicting two-phase  flow  regimes,   experimental  steam-water  data,    and  computer 
sub -routines  for  the  Homogeneous  and  Lockhart-Martinelli  methods  are  given 
in  the  Appendices. 


